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ABSTRACT 


-^7*  Steam-condensation  heat-transfer  measurements  were  made 
using  a  5-tube  in-line  test  condenser  with  an  additional 
perforated  tube  to  simulate  up  to  30  active  tubes.  Results 
were  obtained  for  smooth  tubes,  wire-wrapped  tubes  and 
dropwise-coated  tubes.  The  average  outside  heat-transfer 
coefficient  for  30  smooth  tubes  was  0.64  times  the  Nusselt 
coefficient  for  the  first  tube.  A  total  of  eight  wire- 
diameter  and  wire-pitch  combinations  were  tested:  1.6-mm- 
diameter  wire  wrapped  at  16  mm,  7.6  mm  and  4  mm  wire  pitches, 
1 . 0-mm-diameter  wire  wrapped  at  8  mm,  4  mm  and  2  mm  wire 
pitches,  and  0 . 5-mm-diameter  wire  wrapped  at  4  mm  and  2  mm 
wire  pitches.  The  best  bundle  performance  was  obtained  when 
the  tubes  were  wrapped  with  1. 0-mm-diameter  wire  at  a  wire 
pitch  of  4  mm.  This  combination  resulted  in  an  average  out¬ 
side  heat-transfer  coefficient  for  30  tubes  that  was  1.15 


times  the  value  computed  for  the  first  tube  using  the  Nusselt 


theory.  The  average  outside  heat-transfer  coefficient  for 
the  30  dropwise-coated  tubes  was  1.1  times  the  value  of  the 
heat-transfer  coefficient  for  the  first  tube  in  the  tube 
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bundle.^  Utilizing  either  wire-wrapped  tubes  or  dropwise- 
coated  tubes,  it  is  possible  to  significantly  reduce  the 
condenser  surface  area  and  overall  size. 
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I .  INTRODUCTION 


A.  HISTORICAL  BACKGROUND 

Considerable  interest  has  been  generated  in  reducing  the 
size  and  the  weight  of  propulsion  systems  for  naval  applica¬ 
tions.  Advances  in  condenser  design  could  do  much  to  reduce 
the  size  and  the  weight  of  the  propulsion  plant.  Measures 
to  raise  the  condensing-side  heat-transfer  coefficient  (of 
condenser  tubes)  is  one  way  to  achieve  this  reduction  in 
condenser  size.  This  reduction,  however,  comes  at  a  price. 
This  is  usually  due  to  an  increase  in  the  pumping  power  or 
due  to  an  increase  in  the  initial  cost  of  the  tubes.  For 
naval  applications,  where  the  size  of  a  vessel  may  depend 
upon  the  size  of  the  condenser  (in  a  submarine,  for  example) 
this  reduction  in  the  size  is  justified,  even  at  the  greater 
cost. 

Search  [Ref.  1],  at  the  Naval  Postgraduate  School, 
investigated  the  present  condenser  design  process,  and 
examined  the  potential  benefits  that  might  occur  if  heat- 
transfer  enhancement  was  established  in  the  condenser.  He 
concluded  that  reductions  of  as  much  as  forty  percent  in  the 
size  and  weight  of  condensers  are  possible.  This  is  depen¬ 
dent,  of  course,  on  the  heat-transfer-enhancement  technique 
utilized.  Much  further  research  work  at  the  Naval  Post¬ 
graduate  School  has  been  directed  toward  these  heat-transfer 
enhancement  techniques . 
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Beck  [Ref.  2] ,  Pence  [Ref.  3] ,  Reilly  [Ref.  4] ,  Fenner 
[Ref.  5] ,  and  Ciftci  [Ref.  6]  conducted  research  employing 
a  single-tube  test  condenser.  Their  research  concluded  that 
the  overall  heat-transfer  coefficient  of  enhanced  tubes  may 
be  as  much  as  twice  that  for  smooth  tubes  of  similr  geometry. 
In  a  separate  report,  Marto,  Reilly  and  Fenner  [Ref.  7]  re¬ 
ported  that  most  of  the  increase  in  the  overall  heat-transfer 
coefficient  was  on  the  cooling-water  side  and  was  due  to  an 
increase  of  the  turbulence  and  the  swirl,  as  well  as  to  an 
increase  in  the  inside  surface  area.  Only  a  small  increase 
occurred  on  the  steam  side. 

Present-day  steam  condensers  utilizing  smooth  tubes  are 
limited  in  their  thermal  efficiency,  due  primarily  to  the 
large  thermal  resistances  occurring  on  the  tube  side  of  the 
condenser.  It  is  possible,  however,  by  utilizing  enhanced 
tubes,  to  increase  the  inside  heat-transfer  coefficient  by 
100  percent  or  more.  The  corresponding  increase  in  the  out¬ 
side  heat-transfer  coefficient,  however,  is  less  than  50 
percent.  In  studying  ways  to  further  increase  the  outside 
heat-transfer  coefficient,  Webb  [Ref.  8]  reported  that 
conduction  across  the  condensate  film  is  the  primary  thermal 
resistance  in  film  condensation.  This  thermal  resistance  is 
usually  larger  than  the  thermal  resistance  of  the  tube  wall, 
that  attributed  to  fouling  or  that  due  to  noncondensable 
gases.  It  is  possible  to  reduce  this  thermal  resistance  by 
utilizing  a  geometry  that  reduces  the  film  thickness.  This 


reduction  of  the  thermal  resistance  would  mean  a  correspond¬ 
ing  increase  in  the  outside  heat-transfer  coefficient. 

For  large  tube  bundles,  condensate  inundation  is  present 
and  must  be  considered  when  attempting  to  increase  the 
overall  heat-transfer  coefficient.  Thomas  [Ref.  9]  wrapped 
wire  around  smooth  tubes  in  a  helical  manner  and  tested  the 
condensation  of  ammonia  in  a  large  tube  bundle.  Increases 
of  as  much  as  200  percent  in  the  oustide  heat-transfer 
coefficient  over  that  predicted  by  Nusselt  [Ref.  10]  were 
measured.  This  increase  was  attributed  to  the  effect  of 
surface  tension  drawing  the  condensate  to  the  wire  and 
acting  as  a  condensate  run-off  channel.  In  a  paper  by 
Cunningham  [Ref.  11] ,  "roped"  tubes  were  considered  and 
again  the  increase  in  the  outside  heat-transfer  coefficient 
was  attributed  to  condensate  drainage,  while  an  increase  in 
the  inside  heat-transfer  coefficient  was  attributed  to  the 
increase  in  the  inside  convective  coefficient,  due  to  the 
increased  turbulence.  Kanakis  [Ref.  12]  tested  both  smooth 
and  roped  tubes,  with  and  without  a  wire  wrap  in  an  in-line 
tube  bundle  simulating  up  to  30  tubes.  Adding  the  wire  wrap 
on  the  smooth  tube  increased  the  average,  outside  heat- 
transfer  coefficient  for  the  30  tube  bundle  by  50  percent, 
while  adding  the  wire  warp  to  the  roped  tube  increased  the 
average  outside  heat-transfer  coefficient  for  the  30  tube 
bundle  by  more  than  35  percent. 

For  condenser  tubes,  the  increase  in  the  outside  heat- 
transfer  coefficient  may  also  be  accomplished  by  promoting 


dropwise  condensation.  Tanasawa  {Ref.  13]  reviewed  dropwise 
condensation  and  discussed  the  methods  for  promoting  these 
dropwise  conditions.  He  concluded  that  the  use  of  organic 
polymers  (such  as  Teflon)  was  the  most  promising  of  all  the 
dropwise-promoting  techniques.  Investigations  of  organic 
coatings  by  Brown  [Ref.  14]  have  shown  that  enhancements  of 
up  to  180  percent  are  possible.  Of  primary  concern,  however 
is  that  these  coatings  have  very  low  conductivities  and  must 
therefore  be  ultra-thin.  Another  concern  is  that  these 
coatings  must  be  strongly  adherent  and  sufficiently  tough  to 
withstand  the  conditions  of  their  assembly  and  their  use. 
Holden  [Ref.  15]  investigated  numerous  dropwise-promoting 
coatings.  His  tests  were  based  not  only  on  the  ability  of 
the  coating  to  promote  dropwise  condensation,  but  also  on 
the  ability  of  the  coating  to  sustain  this  dropwise  perfor¬ 
mance  over  an  extended  period  of  time.  In  addition,  the 
coatings  that  were  able  to  sustain  this  dropwise  performance 
were  evaluated  for  their  heat-transfer  performance.  His 
results  indicated  that  the  outside  heat-transfer  coefficient 
for  a  single  tube  could  be  increased  by  a  factor  of  from 
five  to  eight  through  the  use  of  polymer  coatings. 

Since  the  use  of  a  wire  wrap  can  increase  the  heat- 
transfer  performance  of  a  tube  bundle,  an  important  question 
arises:  Is  there  some  optimal  wire  diameter  and  wire  pitch 

combination?  In  addition,  during  dropwise  conditions  what 
is  the  effect  of  condensate  inundation  upon  heat  transfer  in 
a  tube  bundle? 


B.  OBJECTIVES 


i 


i 


The  objectives  of  this  thesis  were  therefore  to: 

1.  Conduct  steam-condensation  tests  to  determine  the 
steam-side  heat-transfer  coefficient  for  16-mm-o.d 
smooth  titanium  tubes, 

2.  Confirm  the  heat-transfer-performance  measurements 
of  Kanakis  [Ref.  12]  on  16-ram-o.d.  smooth  titanium 
tubes  with  a  1.6-mm-o.d.  titanium  wire  wrapped  at 
pitch  of  7.6  mm, 

3.  Conduct  heat-transfer-performance  measurements  to 
determine  the  effect  of  wire  pitch  and  wire 
diameter  on  the  steam-side  heat-transfer 
coefficient  with  inundation,  and 

4.  Conduct  heat-transfer-performance  measurements  to 
determine  the  effect  of  inundation  on  a  dropwise 
coated  tube  and  compare  this  performance  with-  the 
performance  of  a  smooth  tube. 
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II.  THEORETICAL  AND  EMPIRICAL  BACKGROUND 


The  basis  for  the  analysis  of  film  condensation  on  a 
horizontal  tube  was  set  forth  by  Nusselt  in  1916.  His  analy 
sis  was,  however,  for  laminar  film  condensation  on  a  single 
horizontal  tube.  Nusselt' s  analysis  yielded  the  well-known 
relationship  for  the  heat- transfer  coefficient: 


Nu 


=  0.725 


■3  2  . 

kf  Pf  hfQ  g 


Do  ^(Tsat-V 


1/4 


(2 


This  relationship  is  subject  to  the  following  restric¬ 
tions,  as  stated  by  Nobbs  [Ref.  16] : 

1.  The  wall  temperature  is  constant, 

2.  The  flow  is  laminar  in  the  condensate  film, 

3.  The  film  thickness  is  small  compared  to  normal 
tube  diameters, 

4.  All  fluid  properties  are  constant  within  the 
condensate  film, 

5.  Heat  transfer  in  the  film  is  by  conduction, 

6.  All  forces  due  to  hydrostatic  pressure,  surface 
tension,  inertia,  and  vapor/liquid  interfacial  shear 
are  negligible  when  compared  to  viscous  and 
gravitational  forces,  and 

7.  The  vapor/liquid  interface  and  the  surrounding  steam 
are  at  the  saturation  temperature. 

Jakob  [Ref.  17]  extended  the  Nusselt  analysis  to  film 
condensation  on  a  vertical  in-line  column  of  horizontal 
tubes  by  assuming  that  all  the  condensate  from  a  tube  drains 
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Laminar  Flow  Model  Flow  Model 


c.  Eissenberg's  Side-Drainage  Flow  Model 


Figure  2.1  Representations  of  Condensate  Flow 


or  in  terms  of  the  local  heat-transfer  coefficient  for  the 
Nth  tube: 

hN/hx  =  N5/6  -  (N  - 1) 5/6  (2.6 

Eissenberg  [Ref.  19]  did  extensive  experimentation  to 
investigate  the  effects  of  steam  velocity,  condensate  inun¬ 
dation  and  noncondensable  gases  on  the  condensation  heat- 
transfer  coefficient.  He  theorized  that  condensate  does  not 
always  drain  only  onto  tubes  aligned  vertically,  but  can  be 
diverted  sideways,  especially  in  staggered  tube  bundles. 

The  condensate  draining  sideways  strikes  the  lower  tubes  on 
their  sides  rather  than  on  their  tops.  This  is  depicted  in 
Figure  2.1c.  Since  more  heat  is  transferred  from  the  top  of 
a  tube  than  from  its  bottom,  the  net  effect  of  inundation  is 
less  severe.  He  predicted  the  following  relationship: 

h^/h1  =  0.60+0.42  N~1/4  (2.7 

Extensive  experimental  research  into  the  effect  of 
condensate  inundation  has  been  conducted.  However,  the  data 
exhibit  a  substantial  amount  of  scatter  as  shown  in  Figure 
2.2  by  the  cross-hatched  area.  Berman  [Ref.  20]  conducted  a 
compilation  of  film  condensation  data  on  bundles  of  horizontal 
tubes  and  identified  the  following  variables  as  important 
in  causing  the  scatter: 


Figure  2.2  Theories  and  Data  for  Condensate  Inundation 


1.  Bundle  geometry  (in-line  or  staggered) , 

2.  Tube  spacing, 

3.  Type  of  condensing  fluid, 

4.  Operating  pressure, 

5.  Heat  flux,  and 

6.  Local  vapor  velocity. 

Marto  and  Wanniarachchi  [Ref.  21]  noted  that  other  factors 
can  have  an  effect  on  the  data  such  as  noncondensable  gases, 
the  direction  of  the  vapor  flow,  partial  dropwise  condi¬ 
tions,  an  insufficient  abount  of  steam  reaching  the  lower 
tubes  in  the  bundle  and  difficulty  in  the  measurement  of 
the  local  condensing  coefficient. 

Very  small  amounts  of  noncondensable  gases  can  result  in 
significant  reductions  in  the  condensation  heat-transfer 
rate.  Summaries  of  the  phenomenon  have  been  made  by  Chisholm 
[Ref.  22],  and  Webb  and  Wanniarchchi  [Ref.  23],  The  consen¬ 
sus  among  them  is  that  noncondensable  gases  impose  an  added 
thermal  resistance,  since  the  vapor  molecules  must  diffuse 
through  a  gas  layer  prior  to  reaching  the  condensing  surface. 
These  noncondensable  gases  can  also  lead  to  regions  where 
the  tubes  are  inoperative  in  a  condensing  role.  Nonconden¬ 
sable  gases  have  an  adverse  effect  on  the  condenser  perfor¬ 
mance,  and  must  be  taken  into  account  when  designing  a  tube 
bundle . 

The  motion  of  vapor  within  a  condenser  affects  the  film 
condensation  process  because  of  its  effect  on  the  surface 


shear  between  the  vapor  and  the  film,  and  the  resulting  effect 
on  vapor  separation.  Although  the  results  may  be  different 
depending  upon  the  orientation  of  the  steam  flow,  the  general 
effect  of  an  increase  in  vapor  velocity  is  a  corresponding 
increase  in  the  condensing  heat-transfer  coefficient.  There 
exist  a  number  of  both  theoretical  and  empirical  equations 
representing  vapor-shear  effects  on  the  condensate  heat- 
transfer  coefficient  [Refs.  24,  25  and  26]. 

It  is  clear  that  the  measurements  that  were  made  during 
this  thesis  had  both  inundation  and  vapor-shear  effects. 

Marto  [Ref.  27]  stated  that  these  two  effects  are  difficult 
to  separate  from  one  another.  During  this  study,  however, 
an  attempt  was  made  to  separate  these  two  effects.  For  this 
purpose,  data  were  taken  on  the  top  tube  of  the  bundle  with 
vapor  velocity  as  a  variable,  and  the  data  were  correlated 
using  an  expression  similar  to  a  correlation  suggested  by 
Fujii  [Ref.  28] . 

1/2  n 

Nu/Re^  =  BF  (2.8 

This  correlation  is  based  on  numerous  data  for  both 
staggered  and  in-line  tube  bundles.  To  represent  the  test- 
condenser  tube  bundle,  the  constants  B  and  n  in  equation 
(2.8)  were  computed  using  a  least-squares  technique. 


EXPERIMENTAL  APPARATUS 


The  test  facility  was  designed  and  built  originally  by 
Morrison  [Ref.  29] ,  and  modified  by  Noftz  [Ref.  30]  and 
Kanakis  [Ref.  12]  to  simulate  an  active  tube  column  of  up  to 
30  in-line  tubes.  A  detailed  description  of  all  the  com¬ 
ponents  is  contained  in  Ref.  30  and  further  modifications 
are  described  in  Ref.  12.  The  descriptions  in  this  report 
are  brief,  with  particular  focus  given  to  the  modifications 
undertaken  by  the  author. 


A.  STEAM  SUPPLY 

The  steam  supply  system  is  shown  in  Figure  3.1.  House 
steam  flows  through  a  supply  valve  (MS-3) ,  into  a  steam 
separator,  and  is  throttled  down  by  another  valve  (MS-4) 
to  operating  conditions.  The  steam  then  passes  through  an 
orifice  and  into  the  test  condenser  diffuser  before  entering 
the  test  condenser  where  it  flows  through  a  simulated  in-line 
tube  bundle . 


B.  TEST  CONDENSER 

The  dimensions  of  the  test  condenser  shown  in  Figure  3.2 
were  unchanged  from  Noftz ' s  original  design.  The  condenser 
consisted  of  five  active  tubes,  twelve  dummy  tubes  and  one 
perforated  tube.  The  tubes  were  positioned  in  the  test  con¬ 
denser  with  a  pitch-to-diameter  ratio  of  1.5.  The  active 


condensing  length  of  each  tube  was  305  mm.  A  double-walled 
glass  viewing  window  was  provided  on  the  front  of  the  test 
condenser  to  allow  for  visual  observation  of  the  condensation 
process.  Although  the  viewing  window  was  designed  to  have 
heated  air  flowing  between  the  glass  walls,  this  feature  was 
not  utilized  since  no  fogging  of  the  glass  was  observed.  A 
modification  of  the  test  condenser  was  made  to  minimize  the 
steam  flow  along  the  test  condenser  side  walls.  The  modifica¬ 
tion  consisted  of  the  addition  of  twelve  baffle  plates,  placed 
as  illustrated  in  Figure  3.2. 

C.  TEST-CONDENSER  TUBES 

Ten  different  types  of  active  tubes  were  tested  in  this 
set  of  experiments.  The  tubes  were  manufactured  by  Wolverine 
Division  of  Universal  Oil  Products,  and  were  made  of  titanium. 
All  were  of  16  mm  o.d.  and  had  a  minimum  wall  thickness  of 
0.89  mm.  The  ten  types  of  tubes  are  as  listed  below: 

1 .  Smooth  tubes , 

2.  Smooth  tubes  wrapped  with  1.6-mm-o.d.  titanium  wire 
at  three  different  pitches  (16  mm,  7.6  mm  and  4  mm), 

3.  Smooth  tubes  wrapped  with  1.0-mm-o.d.  titanium  wire 
at  three  different  pitches  (8  mm,  6  mm  and  4  mm) , 

4.  Smooth  tubes  wrapped  with  0.5-mm-o.d.  titanium  wire 
at  two  different  pitches  (4  mm  and  2  mm) ,  and 

5.  Smooth  tubes  with  a  dropwise  coating. 

The  actual  wrapping  of  the  wire  around  the  smooth  tubes 
was  performed  by  the  author  of  this  thesis.  A  detailed 
description  of  the  procedure  is  included  in  Chapter  III. 


The  dropwise-coated  tubes  were  sent  to  General  Magna- 
plate  Corporation  where  a  dropwise  coating,  commercially 
referred  to  as  "Nedox,"  was  applied.  First  an  electro- 
deposited  nickel-cobalt  substrate  was  applied  to  the  tube 
surface,  and  then  this  substrate  was  infused  with  a  micro- 
thin  Teflon  layer.  A  heat-treating  cycle  was  employed  to 
ensure  thorough  infusion  of  the  substrate  with  the  Teflon. 

This  coating  is  highly  non-wetting,  has  a  low  coefficient 
of  friction  and  is  heat-resistant.  The  specifics  of  the 
"Nedox"  process  are  proprietory,  but  the  manufacturer  claims 
control  of  the  surface  thickness  to  0.0025  mm.  For-  this  thesis 
the  requested  surface  thickness  was  0.0075  mm. 

D .  PERFORATED  TUBE 

Each  set  of  five  active  tubes  had  its  own  corresponding 
perforated  tube.  This  perforated  tube  was  located  above  the 
uppermost  active  tube;  see  Figure  3.2.  Figure  3.3  shows  a 
typical  pattern  for  the  perforations  of  a  tube,  and  is  an 
illustration  of  the  supply-water  tubing  at  the  test  condenser. 
This  particular  arrangement  of  the  tubing  was  used  in  an 
attempt  to  provide  an  even  flow  from  the  entire  length  of  the 
perforated  tube.  A  heated  tube  was  used  to  provide  the 
water  supply  for  the  perforated  tube.  A  separate  rotameter 
was  provided  in  order  to  control  the  rate  of  flow  of  condensate 
through  the  perforated  tube. 
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Figure  3.3  Plan  View  of  a  Perforated  Tube 


E .  CONDENSATE  SYSTEM 

The  condensate  system  remained  unchanged  from  Kanakis' 
modifications.  Figure  3.1  shows  the  system  as  designed  by 
Noftz  and  as  modified  by  Kanakis.  As  steam  condenses  in  the 
test  condenser,  it  runs  into  and  collects  in  the  test  con¬ 
denser  hotwell.  The  hotwell  is  a  calibrated  cylinder  equipped 
with  a  sight  glass,  providing  a  visual  measure  of  the 
condensate  being  produced.  Measurement  of  the  condensate 
collection  rate  provided  the  flow  rate  to  be  supplied  to 
the  perforated  tube. 

F.  COOLING  WATER  SYSTEM 

Although  the  cooling  water  system  remained  unchanged 
from  Noftz  original  design,  the  five  rotameters  were  recali¬ 
brated  after  Kanakis'  data  were  taken.  Each  active  tube  had 
its  own  rotameter  and  regulating  valve,  and  it  was  possible 
to  control  the  coolant  velocity  through  the  active  tubes  up 
to  values  of  5  m/s.  All  data  for  this  thesis,  however,  were 
taken  at  a  constant  coolant  velocity  of  1.56  m/s.  The  inlet 
temperature  of  the  cooling  water  was  maintained  at  a  nearly 
constant  temperature  by  the  use  of  a  large  supply  tank  and 
a  separate  cooling  tower  to  dissipate  excess  heat  to  the 
atmosphere.  In  order  to  provide  a  bulk  temperature  at  the 
outlet  of  each  active  tube,  a  mixing  chamber  was  installed 
just  upstream  of  the  temperature  monitoring  point. 


G.  INSTRUMENTATION 

1.  Flow  Rates 

Rotameters  were  used  to  measure  the  flow  rate  of  the 
cooling  water  to  each  active  tube  and  to  the  perforated  tube. 

2 .  Temperatures 

Stainless-steel-sheathed  copper-constantan  thermo¬ 
couples  were  used  as  the  primary  temperature  monitoring 
devices.  Table  1  lists  the  locations  monitored  by  each  of 
the  thermocouples.  Calibration,  as  performed  by  Kanakis, 
was  utilized  for  all  data  taken. 

The  perforated-tube  water  supply  was  not  maintainable 
at  a  constant  temperature,  nor  was  it  quickly  adjustable  to 
a  desired  temperature.  The  system,  as  designed,  contained 
a  Gulton  Industries,  West  20  temperature  controller.  Although 
the  manufacturer's  accuracy  was  listed  as  1.25  degrees  F, 
the  system  response  was  so  poor  that  the  temperature  of  the 
perforated-tube  supply  water  varied  by  as  much  as  10  degrees 
C  from  the  beginning  of  the  data  run  to  the  completion  of 
the  data  run. 

3 .  Pressure 

A  Bourdon-tube  pressure  gauge  was  used  to  measure  the 
house  steam  supply  pressure,  while  a  compound  pressure  gauge 
was  used  to  measure  the  steam  pressure  downstream  of  the 
orifice.  A  pressure  transducer  was  used  to  measure  the 
pressure  in  the  test  condenser,  and  a  U-tube  mercury  manom¬ 
eter  was  fitted  to  measure  the  pressure  drop  across  the 


orifice. 


TABLE  1 


Thermocouple  Monitoring  Locations 
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4 .  Data  Collection  and  Display 

A  Hewlett-Packard  model  3054A  Automatic,  Data- 
Acquisition  System,  with  a  HP  2671G  printer  was  used  to 
record  and  display  the  thermocouple  and  transducer  readings. 
The  pressure  transducer  was  assigned  to  channel  19  of  the 
data-acquisition  system,  with  the  thermocouples  assigned  to 
the  channels  as  indicated  in  Table  1. 
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IV.  EXPERIMENTAL  PROCEDURES 


A.  CALIBRATION 

1.  Rotameters 

Calibration  was  performed  on  the  five  rotameters 
that  supplied  the  cooling  water  to  the  active  tubes.  This 
calibration  was  performed  using  a  weighing  tank  and  a  stop¬ 
watch,  and  a  least-squares-f it  calibration  line  was  generated 
for  each  of  the  rotameters. 

2 .  Steam  Orifice 

To  calibrate  the  steam  orifice  the  following  steps 
were  taken: 

1.  Set  a  steam  flow  rate, 

2.  Adjust  the  cooling  water  flow  rate  to  condense  all 
the  incoming  steam, 

3.  Measure  the  pressure  drop  across  the  orifice, 

4.  Measure  the  condensate  collection  rate, 

5.  Repeat  steps  1  through  4  for  different  steam  flow 
rates,  and 

6.  Develop  a  least-squares-f it  straight  line  for  the 
steam  mass  flow  rate  vs.  pressure  drop. 

B.  PREPARATION  OF  CONDENSER  TUBES 

Prior  to  the  winding  of  the  titanium  wire,  the  outside 
surface  of  the  tubes  was  buffed  using  steel  wool  unitl  all 
evidence  of  surface  oxidation  was  removed.  If  any  oxidation 
was  evident  on  the  inner  surface  of  the  tube,  a  test-tube 


brush  was  used  to  apply  a  50-percent  sulfuric-acid  solution 
to  the  inner  surface  of  the  tube.  The  tube  was  then 
thoroughly  rinsed  using  tap  water. 

The  procedure  followed  in  the  wrapping  of  the  wire  onto 
the  smooth  tubes  was  the  same  for  every  wire  diameter  and 
pitch  combination  tested.  The  wire  was  first  welded  to  the 
tube  and  then,  tinder  a  constant  tension  of  22  Newtons,  the 
wire  was  guided  (using  a  metal  plate  to  control  the  spacing 
between  the  wire  wraps)  to  the  proper  pitch  as  the  tube  was 
manually  turned.  Upon  completion  of  the  wrapping,  the  free 
end  of  the  wire  was  then  welded  to  the  tube  to  yield  a 
helically-wrapped  surface  305  mm  in  length.  After  the  wire 
wrapping  of  the  tubes  was  completed  and  just  prior  to  the 
installation  of  the  tubes  in  the  test  condenser,  the  tubes 
were  cleaned  inside  and  out  using  a  brush  and  a  biodegradable 
detergent.  A  thorough  rinsing  was  agin  undertaken,  and  the 
tubes  were  checked  to  ensure  that  they  displayed  proper  wetting 
characteristics.  If  any  non-wetting  areas  were  evident,  a 
50-percent  solution  of  sodium  hydroxide  mixed  with  an  equal 
amount  of  ethyl  alcohol  and  heated  to  about  80  degrees  C,  was 
brushed  onto  the  tube  surface.  After  rinsing  with  tap  water, 
the  tube  was  ready  for  installation  in  the  test  condenser. 

Upon  completion  of  the  installation  of  the  tubes  in  the  test 
condenser,  steam  was  introduced  into  the  test  condenser. 

Again  the  tubes  were  visually  examined.  If  there  were  any 
areas  that  displayed  dropwise  condensation,  the  tubes  were 


r 
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washed  from  above  using  the  perforated- tube  water  supply. 

This  was  sufficient  in  all  cases  to  restore  filmwise 
condensation . 

C.  SYSTEM  OPERATION 

The  system  was  operated  in  accordance  with  the  operating 
instructions  outlined  in  Appendix  A  of  Ref.  17.  The  system 
was  deemed  to  be  operating  at  a  steady-state  condition  when 
the  cooling  water  inlet  temperature  was  steady.  The  cooling 
water  inlet  temperature  was  monitored  using  the  HP-3054A 
Automatic  Data  Acquisition  System,  and  was  considered  steady 
when  two  consecutive  readings  (at  five-minute  intervals) 
varied  by  less  than  0.1  degrees  C. 

One  hour  was  usually  required  from  system  start-up  to 
the  steady-state  condition.  Upon  reaching  the  steady-state 
condition,  a  data  run  was  begun.  Any  changes  in  the  cooling 
water  flow  or  changes  in  the  perforated-tube  water  flow  were 
accompanied  by  a  five-minute  wait  for  system  restabilization, 
prior  to  taking  any  further  data.  The  duration  of  each  data 
set  was  approximately  one  minute ,  and  a  completed  data  run 
consisted  of  30  separate  data  sets.  During  each  data  set, 
the  condensation  collection  rate  was  measured.  For  each 
inundation  condition,  five  consecutive  data  sets  were  taken 
and  the  average  values  were  computed.  After  each  set  of 
five  consecutive  data  sets,  the  average  condensate  collection 
rate  was  used  to  determine  the  rate  of  flow  of  the  water  into 
the  perforated  tube  which  was  to  be  used  for  the  next  set  of 
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data  (the  average  amount  of  condensate  collected  during  a 
data  set  was  introduced  into  the  perforated  tube  for  the 
next  inundation  condition).  For  example,  to  simulate  tubes 
6  through  10,  the  perforated-tube  water  flow  rate  was  set 
equal  to  the  condensate  collection  rate  for  tubes  1  through 
5.  A  data  run  was  considered  completed  upon  the  simulation 
of  inundation  conditions  through  the  30th  tube. 

Thermocouple  readings  and  pressure-transducer  readings 
were  taken  automatically  by  the  data-acquisition  system, 
while  the  settings  of  the  rotameters  and  the  test  condenser 
initial  and  final  hotwell  levels  were  entered  into  the 
computer  using  the  keyboard.  Initial  and  final  hotwell 
levels  were  taken  for  a  time  interval  of  one  minute. 

The  following  conditions  were  the  operating  conditions 
for  all  data  taken  for  this  thesis: 

1.  The  coolant  velocity  was  1.55  ±  0.05  m/s, 

2.  The  coolant  inlet  temperature  was  24  degrees  C.  (±3 
degrees  C.  depending  on  the  day  of  the  run,  but 
constant  during  a  given  run) ,  and 

3.  The  saturation  temperature  was  100.5  ±  0.5  degrees  C. 
D.  DATA- REDUCTION  PROGRAM 

A  computer  program  was  utilized  to  process  and  plot  all 
raw  data.  The  program  was  written  in  BASIC  language  and 
was  run  on  an  HP-9826  computer.  A  listing  of  this  program 
is  included  as  Appendix  C. 


I 


38 


-  ■-  I  ■  J  ",  \  •  » 


E.  HEAT-TRANSFER-COEFFICIENT  CALCULATION 

1.  Inside  Heat-Transfer  Coefficient 
The  inside  heat-transfer  coefficient  was  computed 

using  the  Sieder-Tate  equation  described  in  Holman  [Ref. 

31]  : 

Nu  =  h±  D./kc  =  Ci  Re0*8  Pr1/3  (i^/pj0*14  (4.1) 

where  the  coefficient  C^  was  computed  using  a  modified 
Wilson  plot  [Ref.  12] .  Using  this  technique,  Kanakis  [Ref. 

12]  obtained  a  Sieder-Tate  coefficient  of  0.029  ±  0.001. 

The  method  Kanakis  used,  however,  did  not  include  the  variation 
of  the  condensate  film  properties  as  a  function  of  heat  flux. 

A  slightly  different,  modified  Wilson  plot,  as  described  by 
Nobbs  [Ref.  16] ,  that  does  account  for  this  variation  in 
the  film  properties,  was  used  in  this  thesis.  Appendix  D 
gives  a  description  of  this  method.  Reprocessing  Kanakis* 
data  with  this  method,  the  Sieder-Tate  coefficient  was  deter¬ 
mined  to  be  0.028  ±  0.001. 

2 .  Outside  Heat-Transfer  Coefficient 
The  heat-transfer  rate  to  the  cooling  water  can  be 

computed  using  an  energy  balance,  and  can  be  expressed  in 
terms  of  an  overall  heat-transfer  coefficient  by: 


Q 


m 


pc 


(T 


CO 


-T  .  ) 
ci 


U  A  LMTD 
o  o 


(4.2) 


where 


LMTD 


(4.2a) 


T  -T  . 
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Solving  for  the  overall  heat-transfer  coefficient  gives: 
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=  m 


c  /A 
pc'  o 


In  ( (T 


sat 
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-Tco>> 


(4.3) 


The  overall  thermal  resistance  can  be  written  as  the  sum  of 
the  internal,  wall  and  external  resistances: 


U  A 
o  o 


i  +  j*  +  _  L_ 

h.  A.  A  h  A 
11  o  o  o 


where 


(4.4) 


R»  -  Do  ln<VDi1''2  km 


(4 .4a) 


or  from  this  result,  the  outside  heat-transfer  coefficient 
may  be  written: 


1/U -  D  /D.  h.  -  R 
'  O  o'  1  1  w 


(4.4b) 


F.  DATA-REDUCT ION  TECHNIQUE  FOR  THE  OUTSIDE  HEAT-TRANSFER 
COEFFICIENT 

A  listing  of  the  computer  program  used  in  the  calculation 
of  the  outside  heat-transfer  coefficient  can  be  found  in 
Appendix  C.  The  steps  required  in  the  technique  are  outlined 
below: 
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Calculate  the  average  bulk  cooling  water  temperature 


T.  =  +  T  . )/2 

b  co  ci 


Calculate  the  cooling  water  velocity! 


V  =  m  /A.  p„ 
C  C  1  c 


Calculate  the  cooling  water  Reynolds  number; 


“  ”  “c  Vc  °i/vc 


Calculate  the  heat  transferred  to  the  cooling 
water  using  the  left-hand  side  of  equation  (4.2) 


Calculate  the  heat  flux: 


q  =  Q/(tt  Dq  L) 


Since  the  film  temperature  was  not  a  known  quantity, 
an  iterative  scheme  was  employed  to  calcualte  the 
film  temperature,  as  indicated  below: 


a) 


Assume  a  film  temperature 


(say  Tf  = 


Tsat> 


9 


b)  Calculate  the  Nusselt  coefficient. 


Nu 


=  0.651 


-.3  2  . 

kf  Pf  hfq  9 

yf  Do  q 


1/3 


c)  Evaluate  T,  ,  and 


d)  Repeat  steps  a)  through  c)  until  T.  is 
approximately  equal  to  T^.  t,c 

Calculate  the  inside  heat- transfer  coefficient: 


hi  *  kc  Ci  Re°'8  prl/3  Cf/Di 


where  is  the  correction  factor  (v>c/y  ) 


0.14 


(4.11) 


(4.11a) 


Since  the  correction  factor  is  evaluated  at  the 
inner  wall  temperature,  an  iterative  scheme  was 
employed  to  calculate  the  wall  temperature  as 
indicated  below: 

a)  Assume  a  value  for  Cf  (e.g.,  1.0), 

b)  Calculate  h^  using  equation  (4.11), 

c)  Calculate  a  water-side  temperature  drop. 


DT  =  q  D./h.  Dq 


d)  Calculate  a  new  correction  factor,  and 


(4.12) 


C-.  =  (u  /u  ) 

f,c  VMc/Mw 


0.14 


(4.13) 


e)  Repeat  steps  a)  through  d)  until  Cf  is 
approximately  equal  to  Cj.  ' 

Calculate  the  log-mean-temperature  difference  using 
equation  (4.2a). 

Calculate  the  overall  heat-transfer  coefficient: 


11.  Calculate  the  outside  heat-transfer  coefficient  using 
equation  (4.4b). 

12.  Calculate  the  normalized,  local,  outside  heat- 
transfer  coefficient. 


13.  Calculate  the  normalized,  average,  outside  heat- 
transfer  coefficient. 

vhi 

NOTE:  In  the  above  formulation, the  thermal  resistances  due 

to  noncondensable  gases  and  any  fouling  were  assumed  to  be 
negligible . 


G.  DATA- REDUCTION  TECHNIQUE  FOR  THE  VAPOR-SHEAR  CORRELATION 


For  the  calculation  of  the  vapor-shear  correlation, 
measurements  on  only  the  first  active  tube  in  the  tube 
bundle  were  taken. 

1.  Calculate  the  inside  heat-transfer  coefficient  as 
outlined  above. 

2.  Calculate  the  outside  heat-transfer  coefficient 
as  outlined  above. 

3.  Calculate  a  temperature  correction: 


DT  =  q/hQ 


(4.15) 


4.  Assume  a  film  temperature: 


T 


T 


DT/2 


(4.16) 


Calculate  the  steam  velocity: 


Vv  =  mv  VAf 


where 

Aj.  =  2  L(Pt  Pl-»  D^/4)/Pl 


Af  is  the 
[Ref.  16] 

Calculate 


Mean  Vapor  Flow  Area  defined  by  Nobbs 
and  depicted  in  Figure  4.1. 

the  two-phase  Reynolds  number: 


Re24,  =  pf  Vv  Do/vf 


Calculate  the  dimensionless  quantity  F: 


F  *  Do  vf  hfg  kf  DT 


Calculate  the  Nusselt  number: 


Nu  =  hQ  Do/kf 


(4.17) 


(4.17a) 


(4.18) 


(4.19) 


(4.20) 


Calculate  Y: 


To  calculate  the  mean  flow  area  the  following  steps 
are  required: 

1)  Calculate  the  cross-hatched  area  above, 
AREA  =  PT  PL  -  IT  Dq2  /  4 

2)  Divide  by  longitudinal  pitch  to  obtain  an 
average  flow  dimension 

Dp  =  AREA  /  PL 

3)  Multiply  by  the  condenser  length  and  the 
number  of  flow  paths  (2  flow  paths  in  the 
case  of  this  thesis). 

Ap  -  2  L  (PT  PL  -  T  D02  /  4)  /  PL 

Figure  4.1  Mean  Vapor  Flow  Area 
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0. 


Repeat  steps  1  through  9  at  different  steam  mass 
flow  rates  into  the  test  condenser.  For  this  thesis 
10  different  steam  mass  flow  rates  were  considered, 
and  three  measurements  taken  at  each  mass  flow  rate. 


L.  Starting  with  an  assumption  that  the  data  would 
be  of  a  form  similar  to  that  postulated  by  Fujii 
[Ref.  28] ,  an  expression  was  assumed: 


Y  =  B  Fn;  F  as  defined  in  equation  (4.19)  (4.22) 


2.  Perform  a  least-squares-fit  of  the  data  and 
determine  the  constant  B  and  the  exponent  n. 
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RESULTS  AND  DISCUSSION 


A.  RESULTS  BEFORE  AND  AFTER  TEST  CONDENSER  MODIFICATION 

Initially,  data  were  taken  with  the  test  condenser  used 
by  Kanakis  [Ref.  12] ,  with  no  modifications.  As  can  be  seen 
from  Figure  5.1,  the  nondimens ionalized  heat-transfer  coeffi¬ 
cient  exhibits  a  saw-toothed  variation  over  the  inundation 
range  considered.  The  particular  wire  pitch  and  wire  diameter 
combination  chosen  for  Figure  5.1  was  representative  of  the 
pattern  apparent  in  the  data  sets  taken  prior  to  the  test- 
condenser  modification. 

It  was  determined  that  the  original  design  of  the  test 
condenser  allowed  too  much  steam  to  bypass  the  five  active 
tubes,  thereby  allowing  the  steam  to  flow  between  the  walls 
of  the  test  condenser  and  the  two  rows  of  dummy  tubes.  Since 
steam  was  bypassing  the  active  tubes,  inadequate  steam  was 
flowing  in  the  vicinity  of  the  active  tubes.  This  led  to 
artificially  low  heat-transfer  readings  in  the  lower  tubes 
of  the  tube  bundle.  To  minimize  this  problem,  baffle  plates 
were  installed  on  each  side  of  the  test  condenser  as  des¬ 
cribed  previously  (Figure  3.2). 

After  the  placement  of  the  baffle  plates,  the  previously 
tested  tube  sets  were  retested.  As  can  be  seen  in  Figure 
5.1,  the  saw-toothed  variation  of  the  curve  is  much  less 
pronounced  and  appears  to  be  nearly  corrected.  The  net 
effect  of  the  removal  of  the  saw-toothed  variation  by  the 


Figure  5.1  Data  With  and  Without  Baffles 


installation  of  the  baffle  plates  is  an  increase  in  the 
measured  heat-transfer  coefficient.  This,  in  turn,  means 
that  the  effect  of  inundation  is  not  actually  as  great  as 
had  been  earlier  suspected.  All  data  represented  in  this 
report  will,  hereafter,  be  with  the  baffles  in  place. 

B.  VAPOR-SHEAR  CORRELATION 

As  discussed  previously,  the  vapor  shear  correlation  was 
assumed  to  be  of  the  form  postulated  by  Fujii  [Ref.  28] — 
see  equation  (4.8).  Figure  5.2  is  the  graph  of  Y  vs.  F  for 
the  first  active  smooth  titanium  tube.  For  this  tube,  the 
value  of  B  was  determined  to  be  0.834  and  the  value  of  n  was 
determined  to  be  0.166.  These  results  fall  about  15  percent 
below  the  Fujii  correlation  [Ref.  28]  (B  =  0.96  and  n  =  0.20) 


There  are  several  possible  explanations  for  this  difference. 


1.  Both  the  value  of  Y  and  the  value  of  F  are  partly 
based  on  the  steam  velocity,  with  a  higher  value  of 
steam  velocity  resulting  in  a  lower  value  for  both  Y 
and  F.  During  this  thesis,  the  steam  velocity  was 
only  varied  from  1.4  to  2.1  m/s.  The  corresponding 
range  of  F  was  therefore  from  3.6  to  1.5,  while  the 
corresponding  range  of  Y  was  from  1.1  to  0.8.  Since 
the  range  of  both  F  and  Y  were  considerably  more  exten¬ 
sive  for  the  calculation  of  the  Fujii  correlation, 
this  is  a  possible  explanation  for  the  discrepancy. 

2.  For  the  Fujii  correlation,  both  in-line  and  staggered 
tube  bundles  were  considered.  This  could  also  be  a 
possible  reason  for  the  discrepancy. 

3.  The  direction  of  the  steam  flow  considered  for  the 
determination  of  the  Fujii  correlation  included 
vertically  downward  flow,  vertically  upward  flow  and 
also  horizontal  flow.  This  difference  in  flow 
directions  is  another  possible  explanation  for  the 
discrepancy  between  the  correlation  determined  in 
this  report  and  the  Fujii  correlation. 
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Figure  5.2  Vapor-Shear  Correlation 


4.  As  mentioned  previously,  baffle  plates  were  installed 
in  order  to  reduce  the  steam  bypassing  the  active 
tubes.  Even  though  the  baffles  reduced  the  amount  of 
steam  bypassing  the  active  tubes,  this  is  no  guarantee 
that  the  steam  still  able  to  bypass  the  active  tubes 
could  not  have  had  some  influence  on  the  vapor-shear 
correlation  calculated  in  this  thesis. 

Since  Reference  28  does  not  specifically  say  how  the 
flow  area  used  in  his  correlation  was  determined,  this  is 
another  possible  area  for  discrepancy.  In  addition,  Fujii 
[Ref.  28]  does  not  show  which  tube  in  the  bundle  was  used 
for  the  formulation  of  his  correlation.  The  actual  location 
would  have  an  important  effect  on  the  vapor  velocity  and 
could  be  another  possible  area  for  discrepancies  to  develop. 

Vapor-shear  correlations  were  determined  for  all  of  the 
tube  sets  tested  for  this  report.  Table  2  is  a  summary  of 
the  constants  determined  for  the  vapor-shear  correlations. 

C.  THE  EFFECT  OF  INUNDATION  ON  A  SMOOTH  TUBE  BUNDLE 

1 .  With  Vapor-Shear  Contribution  Present 

Figure  5.3  shows  the  variation  of  the  normalized, 
local  heat-transfer  coefficient  for  up  to  30  tubes.  The 
value  for  the  normalized,  local  heat-transfer  coefficient 
for  the  30th  tube  is  approximately  45%  above  the  value  pre¬ 
dicted  by  Nusselt  theory,  but  very  near  the  value  predicted 
by  the  Kern  relationship.  As  discussed  in  Chapter  II,  the 
Nusselt  theory  for  a  tube  bundle  is  based  on  idealized  laminar 
flow  from  tubes  above  to  tubes  below.  The  idealizations  of 
Nusselt  result  in  a  conservative  estimate  of  the  heat-transfer 
coefficient,  so  the  fact  that  the  data  are  45%  above  the 


TABLE  2 


Vapor-Shear  Correlations 


TUBE  WIRE 


SET 

DIAMETER  (mm) 

PITCH  (mm) 

_B 

n 

1 

Smooth  tube 

0.834 

0.185 

2a 

1.58 

16 

0.835 

0.178 

b 

1.58 

7.6 

0.837 

0.189 

c 

1.58 

4 

0.725 

0.192 

3a 

1.01 

8 

0.916 

0.161 

b 

1.01 

6 

0.871 

0.173 

c 

1.01 

4 

0.876 

0.220 

4a 

0.50 

4 

1.036 

0.200 

b 

0.50 

2 

0.915 

0.192 

5  . 

Dropwise-coated  tube 

0.796 

0.166 
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Figure  5.3  Inundation  Effect  on  a  Smooth  Tube  Bundle 

(Normalized,  Local  Heat-Transfer  Coefficient) 


Nusselt  prediction  is  not  totally  unexpected.  Since  the  Kern 
relationship  is  a  more  realistic  view  of  the  condensate 
action  in  the  tube  bundle,  one  would  expect  the  data  to  be 
much  nearer  to  the  Kern  relationship — as  it  is.  The  close 
agreement  with  the  Kern  relationship  also  was  considered  to 
be  an  indication  that  the  test  apparatus  were  operating 
properly. 

Figure  5.4  shows  the  normalized,  average  heat-transfer 
coefficient  for  up  to  30  tubes.  The  curve  generated  by  these 
data  is  much  smoother  than  the  curve  generated  by  the 
normalized,  local  heat-transfer  coefficients.  The  value  for 
the  normalized,  average  heat-transfer  coefficient  of  the  40th 
tube  lies  about  26%  above  the  value  predicted  by  Nusselt, 
but  lies  only  5%  below  the  value  predicted  by  Kern.  Assuming 
a  generalized  form  of  the  average  heat-transfer  coefficient 
for  N  tubes,  divided  by  the  heat-transfer  coefficient  for 
the  first  tube  in  the  bundle,  yields: 

h^/h1  =  N~S  (5.1 


Referring  to  Figure  5.4,  the  dashed  line  roughly  following 
the  data  represents  a  least-squares-f it  exponential  curve 
for  the  data,  and  the  exponent  derived  is  s  =  0.183.  This 
exponent  and  the  exponent  suggested  by  Kern  (s  =  0.167)  are 
in  reasonable  agreement. 
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Figure  5.4  Inundation  Effect  on  a  Smooth  Tube  Bundle 

(Normalized,  Average  Heat-Transfer  Coefficient) 


.  With  Vapor-Shear  Contribution  Removed 


Figure  5.3  shows  the  variation  of  the  normalized, 
local  heat-transfer  coefficient  for  up  to  30  tubes  with  the 
contribution  due  to  vapor  shear  removed.  The  vapor-shear 
correlation  used  was  the  one  presented  earlier  in  this  report 
(B  =  0.834  and  n  =  0.185).  As  can  be  seen  from  the  graph, 
the  trend  in  the  data  has  now  become  a  reversed  saw-tooth. 
This  is  an  indication  of  overcorrecting  for  the  effect  of 
vapor  shear.  This  unacceptable  trend  in  the  data  shows  that 
mathematical  separation  of  the  vapor-shear  contribution  is 
not  possible  as  stated  by  Marto  [Ref.  27] .  Therefore, 
throughout  the  remainder  of  this  thesis  the  vapor-shear  con¬ 
tribution  is  included,  and  no  further  attempt  is  made  to 
remove  the  contribution  due  to  vapor  shear. 

D.  THE  EFFECT  OF  INUNDATION  ON  SMOOTH  TUBES  WRAPPED  WITH 

1. 6-MM-O.D.  WIRE 

Even  though  three  different  pitches  (16  mm,  7.6  mm  and 
4  mm)  were  tested  at  this  wire  diameter,  only  the  graph  for 
the  optimum  pitch  will  be  presented  for  the  normalized, 
local  heat-trar.sfer  coefficient  for  a  bundle  of  up  to  30 
tubes.  The  graph  for  the  pitch  of  7.6  mm  is  presented  in 
Figure  5.5.  The  value  for  the  normalized,  local  heat- 
transfer  coefficient  for  the  30th  tube  is  approximately  170% 
above  the  value  predicted  by  the  Nusselt  theory  and  approxi¬ 
mately  83%  above  the  value  predicted  by  the  Kern 
relationship. 


Figure  5.5  Effect  of  1.6-mm-o.d.  Wire  on  Inundation 

(Normalized,  Local  Heat-Transfer  Coefficient) 


Figure  5.6  shows  the  normalized,  average  heat-transfer 
coefficient  for  a  bundle  of  up  to  30  tubes.  All  three 
pitches  are  represented,  with  the  values  of  the  normalized, 
average  heat-transfer  coefficient  of  the  30th  tube  for  the 
pitches  of  16,  7.6  and  4  mm  approximately  45%,  118%,  and 
120%,  respectively,  above  the  value  predicted  by  Nusselt. 

As  mentioned  previously,  the  wire  wrapping  makes  two  separate 
contributions  towards  the  increase  in  the  heat-transfer 
coefficient — it  thins  the  condensate  film  and  acts  as  a 
condensate  drainage  channel.  When  compared  with  Figure  5.3, 
Figure  5.5  shows  that  condensate  inundation  has  less  effect 
on  wire-wrapped  tubes  than  it  does  on  a  set  of  smooth  tubes 
and  that  a  wire  pitch  of  7.6  mm  gives  the  optimal  results. 

The  optimum  wire  pitch  or  wire  diameter  is  not  simply 
based  on  the  largest  value  of  but  also  on  the  actual 

value  of  h^.  For  example,  both  4  mm  wire  pitch  and  7.6  mm 
wire  pitch  resulted  in  about  the  same  1^/1^  value  (see 
Figure  5.6).  But  the  h^  on  the  tube  wrapped  with  the  7.6  mm 
wire  pitch  is  about  10  percent  greater  than  the  tube  with 
the  4.0  mm  wire  pitch.  Thus,  the  7.6  mm  wire  pitch  gave 
the  optimal  results  for  the  1.6-mm-o.d.  wire.  A  more 
thorough  discussion  of  the  point  is  provided  later  in  this 
chapter  (Section  H) . 

E.  THE  EFFECT  OF  INUNDATION  ON  SMOOTH  TUBES  WRAPPED  WITH 

1.0-MM-O.D.  WIRE 

Three  different  pitches  (8  mm,  6  mm  and  4  mm)  were 
tested  at  this  wire  diameter.  Again,  only  the  graph  for  the 
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optimum  pitch  will  be  presented  for  the  normalized,  local 
heat-transfer  coefficient  for  a  bundle  of  up  to  30  tubes. 

The  graph  for  a  pitch  of  4  mm  is  presented  as  Figure  5.7. 

The  value  for  the  normalized,  local  heat-transfer  coefficient 
for  the  30th  tube  is  approximately  151%  above  the  value 
predicted  by  the  Nusselt  theory  and  71%  above  the  value 
predicted  by  the  Kern  relationship. 

Figure  5.8  shows  the  normalized,  average  heat-transfer 
coefficient  for  a  bundle  of  up  to  30  tubes.  Again,  all 
three  pitches  are  represented.  The  values  of  the  normalized, 
average  heat-transfer  coefficient  of  the  30th  tube  (for  the 
pitches  8,  6  and  4  mm)  are  approximately  90%,  105%  and  117%, 
respectively,  above  the  value  predicted  by  Nusselt.  These 
graphs  again  show  that  this  tube  configuration  is  less 
affected  by  condensate  inundation  than  the  smooth-tube  set, 
and  that  the  optimal  wire  pitch  is  4  mm  with  this  smaller 
diameter  wire. 

F.  THE  EFFECT  OF  INUNDATION  ON  SMOOTH  TUBES  WRAPPED  WITH 

0 . 5-MM-O.D.  WIRE 

Two  pitches  (4  mm  and  2  mm)  were  tested  at  this  wire 
diameter.  For  the  normalized,  local  heat-transfer  coeffi¬ 
cient  for  a  bundle  of  up  to  30  tubes,  only  the  optimum  pitch 
of  4  mm  will  be  presented.  The  graph  is  presented  as  Figure 
5.9.  The  value  for  the  normalized,  local  heat-transfer 
coefficient  for  the  30th  tube  is  approximately  101%  above 
the  value  predicted  by  Nusselt  and  37%  above  the  value 
predicted  by  the  Kern  relationship. 


Figure  5.7  Effect  of  1.0-mm-o.d.  Wire  on  Inundation 

(Normalized,  Local  Heat-Transfer  Coefficient) 


Figure  5.8  Effect  of  1.0-inm-o.d.  Wire  on  Inundation 

(Normalized,  Average  Heat-Transfer  Coefficient) 


Figure  5.9  Effect  of  0 . 5-mm-o .d .  Wire  on  Inundation 

(Normalized,  Local  Heat-Transfer  Coefficient) 


Figure  5.10  shows  the  normalized,  average  heat-transfer 
coefficient  for  a  bundle  of  up  to  30  tubes.  Both  of  the 
pitches  are  represented.  The  values  of  the  normalized, 
average  heat-transfer  coefficient  for  the  30th  tube  (for  the 
pitches  4  and  2  mm)  are  approximately  71%  and  80%,  respectively 
above  the  value  predicted  by  Nusselt.  These  graphs  also 
indicate  that  this  tube  configuration  is  not  affected  by 
condensate  inundation  as  much  as  the  smooth-tube  set. 

G.  THE  EFFECT  OF  INUNDATION  ON  SMOOTH  TUBES  WITH  A  DROPWISE 

COATING 

A  set  of  smooth  tubes  with  a  dropwise  coating  was  tested 
primarily  to  observe  the  effect  of  inundation.  Figure  5.11 
shows  the  normalized,  average  heat-transfer  coefficient  for 
a  bundle  of  30  tubes.  Also  included  in  the  figure  for 
comparison  are  the  data  for  the  smooth  tube.  As  is  immedi¬ 
ately  evident  from  the  graph,  condensate  inundation  does  not 
decrease  the  performance  of  the  dropwise-coated  tubes.  The 
ratio  of  the  average  outside  heat-transfer  coefficient  of  N 
tubes,  divided  by  the  heat- transfer  coefficient  of  the  first 
tube  can  be  seen  to  rise  at  first  and  then  remain  steady  at 
a  value  of  approximately  1.1. 

Since  the  coating  is  non-wetting,  no  large  droplets  form; 
rather,  the  droplets  are  small  when  they  first  form  and  roll 
down  the  side  of  the  tube,  as  soon  as  they  reach  a  moderate 
size.  As  these  droplets  are  making  their  descent  down  the 
side  of  the  tube,  they  collect  additional  droplets  until  they 


(Normalized,  Average  Heat-Transfer  Coefficient) 


reach  a  point  where  they  are  able  to  fall  off  the  tube.  As 
these  droplets  strike  the  tubes  below,  they  again  collect 
condensate  and  act  as  a  sweeper  of  the  condensate  across  the 
surface  of  the  tubes  below. 

H .  SUMMARY 

Table  3  represents  a  summary  of  the  results  obtained 
from  the  data  sets  presented  in  this  thesis.  The  column 
hi/hNu  is  the  measure  of  the  local  heat-transfer  coefficient 
of  the  first  tube  compared  to  the  Nusselt  prediction.  It 
can  be  considered  a  measure  of  the  effectiveness  of  the 
wire-and-pitch  combination  to  thin  the  condensate  film 
between  successive  wire  wraps.  Refer  to  Figure  5.12  for 
clarification  of  this  point.  The  column  hj^/h^  is  the 
measure  of  the  average  heat-transfer  coefficient  for  30 
tubes  compared  to  the  coefficient  for  the  first  tube.  This 
can  be  considered  to  be  a  measure  of  the  effectiveness  of 
the  wire-and-pitch  combination  at  handling  the  condensate 
drainage.  The  column  ^q/^nu  t*le  measure  of  the  average 
heat-transfer  coefficient  for  30  tubes  compared  to  the  heat- 
transfer  coefficient  predicted  by  the  Nusselt  theory  for  a 
single  tube.  It  can  be  viewed  as  a  measure  of  the  overall 
performance  of  the  30  tube  bundle.  The  last  column  (labelled 
s)  is  the  exponent  defined  in  equation  (5.1),  and  is  calcu¬ 
lated  using  the  least-squares  fit  for  each  data  run. 

As  can  be  seen  from  the  data,  it  appears  that  the  larger 
the  wire,  the  better  it  aids  in  condensate  drainage. 
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Data  Summary 

TUBE  WIRE  WIRE  2 

SET  DIAMETER  (nun)  PITCH  (mm)  h,  (kW/m  .  K)  l^/h 


NOTE:  Determined  by  curve  fitting  the  data  for  tubes  11  through  30. 


Conversely,  the  smaller  the  wire  the  higher  the  heat-transfer 
performance  of  the  first  tube.  These  last  two  observations 
are  based  on  a  comparison  at  a  4  mm  pitch  but  with  three 
different  wire  diameters.  The  explanation  for  this  is  that, 
for  a  given  pitch,  as  the  wire  diameter  of  the  wrap  is  in¬ 
creased,  less  and  less  tube  surface  area  remains  exposed  to 
the  steam  flow.  In  addition  to  this,  the  condensate  film 
thickness  may  also  increase  as  the  wire  diameter  is  increased. 
Thus,  for  a  given  wire  pitch,  the  heat-transfer  coefficient 
of  the  first  tube  will  increase  with  a  decrease  in  the  wire 
diameter.  On  the  other  hand,  the  amount  of  condensate  col¬ 
lected  around  the  wire  increases  with  increasing  wire  diameter 
(see  Figure  5.12) .  While  the  surface-tension  forces  tend  to 
retain  liquid  at  the  wire,  gravitational  forces  tend  to  draw 
liquid  along  the  wire.  As  the  wire  diameter  increases,  the 
gravitational  forces  increase  at  a  faster  rate  than  the 
surface-tension  forces.  Therefore,  condensate  drainage  in¬ 
creases  with  increasing  wire  diameter. 

I .  OBSERVATIONS 

1 .  Smooth  Tubes 

During  the  data  runs,  filmwise  condensation  was 
observed  with  no  visible  indication  of  dropwise  condensation. 
When  the  inundation  rate  was  increased,  more  drops  formed 
along  the  bottom  of  the  tubes,  but  the  droplet  size  did  not 
increase  noticeably.  However,  at  the  higher  inundation 


As  the  steam  was  condensing  on  the  tubes,  the  conden¬ 
sate  was  observed  to  be  drawn  to  the  wire  wrap,  leaving  a 
thin  film  between  the  wires  (film  thickness  was  dependent 
upon  the  particular  wire  diameter  and  pitch) .  The  condensate 
then  ran  down  the  wire  and  formed  droplets  on  the  bottom  of 
the  tube.  The  droplet  size  was  different  for  each  configura¬ 
tion  of  wire  diameter  and  pitch.  For  the  larger  diameter 
wire,  rivulets  formed  and  there  was  a  solid  condensate 
"bridge"  from  the  tubes  above  to  the  tubes  just  below  them. 

As  the  inundation  rate  increased,  these  bridges  became  more 
pronounced.  It  was  also  noticed  that  the  tendency  to  splash 
increased  as  the  inundation  rate  increased,  and  was  quite 
noticeable  for  all  of  the  tube  sets  tested. 

3 .  Droowise-Coated  Tubes 


As  was  mentioned  previously,  the  condensate  droplets 
from  the  upper  tubes  in  the  bundle  act  as  sweepers  of  the 
condensate  across  the  surface  of  the  tubes  below  them  in  the 
tube  bundle.  During  the  data  run,  the  condensate  collecting 
on  the  upper  tubes  and  then  falling  onto  the  tubes  below  was 
sweeping  the  lower  tubes  at  such  a  rate  that  the  droplets 
forming  on  the  lower  tubes  never  reached  a  size  larger  than 
approximately  2  mm. 
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VI.  CONCLUSIONS 


The  average,  outside  heat-transfer  coefficient  for 
30  smooth  tubes  in  a  vertical  column  was  64%  of  the 
Nusselt  coefficient  for  the  first  tube  in  the  tube 
bundle. 

With  wire  wrapping,  the  average,  outside  heat-transfer 
performance  of  the  tube  bundle  is  considerably  improved 
Significant  improvements  in  the  outside  heat-transfer 
performance  are  also  possible  by  utilizing  dropwise 
coatings . 

The  optimum  pitch  for  the  1.6-mm-o.d.  wire  wrapping 
occurred  at  7.6  mm.  With  this  wire  arrangement,  the 
average,  outside  heat-transfer  coefficient  for  the 
30  tubes  was  92%  of  the  Nusselt  coefficient  for  the 
first  tube  in  the  tube  bundle. 

The  best  pitch  tested  for  the  1.0-inm-o.d.  wire 
wrapping  occurred  at  4  mm.  With  this  wire  arrangement, 
the  average,  outside  heat-transfer  coefficient  for 
the  30  tubes  was  115%  of  the  Nusselt  coefficient  for 
the  first  tube  in  the  tube  bundle. 

The  best  pitch  tested  for  the  0.5-mm-o.d.  wire  wrapping 
occurred  at  4  mm.  In  this  case,  the  average,  outside 
heat-transfer  coefficient  for  the  30  tubes  was  110% 
of  the  Nusselt  coefficient  for  the  first  tube  in 
the  tube  bundle. 

In  simulating  a  bundle  of  30  tubes  coated  with  a 
dropwise  promoter,  the  normalized,  average,  outside 
heat-transfer  coefficient  increased  by  approximately 
10%  over  the  value  for  the  top  tube  in  the  bundle. 


VII.  RECOMMENDAT IONS 


A.  TEST  APPARATUS  MODIFICATIONS 

The  modifications  listed  below  should  be  considered 
prior  to  continued  use  of  the  test  apparatus: 

1.  Redesign  the  test  condenser  to  reduce  the  possibility 
of  the  steam  bypassing  the  active  test  condenser  tubes 

2.  Install  a  larger  heating  system  in  the  perforated- 
tube  water  supply  tank  or  modify  the  existing  heating 
system  to  increase  the  slow  response  of  the  heating 
system. 

The  modifications  listed  below  would  be  beneficial  to 
the  operation  of  the  test  apparatus  but  are  not  necessary 
for  proper  operation: 

1.  Redesign  the  existing  controls  so  that  the  test 
apparatus  is  controllable  from  a  single  panel. 

2.  Redesign  the  test  condenser  hotwell  to  provide  a  more 
accurate  measure  of  the  condensate  collection  rate. 


B.  ADDITIONAL  TESTS  TO  CONDUCT 

The  following  additional  tests  would  be  important  in 
the  continuation  of  this  investigation: 

1.  Conduct  tests  with  additional  wire  pitches  for  the 
1 . O-mm-diameter  wire  to  determine  the  optimum  wire 
pitch. 

2.  Conduct  tests  using  other  dropwise  coatings. 

3.  Conduct  tests  with  the  commercially-available  finned 
tubes . 

4.  Conduct  additional  tests  to  determine  the  effect  of 
vapor  shear  on  the  outside  heat-transfer  coefficient. 


APPENDIX  A 


SAMPLE  CALCULATIONS 


A.  A  sample  calculation  is  performed  in  this  section  to 
illustrate  the  solution  procedure  used  in  the  data-reduction 
program  presented  in  Appendix  C.  The  calculations  are  limited 
to  the  first  tube  in  the  bundle  only.  All  thermophysical 
properties  were  determined  from  the  Tables  in  Reference  32. 


EXPERIMENTAL  CONDITIONS: 

Pressure  Condition 

Inundation  Condition 

Inlet  Temperature  of  Cooling  Water 

Outlet  Temperature  of  Cooling  Water 

Saturation  Temperature 

Cooling  Water  Rotameter  Setting 


Atmospheric 
5  tubes 
24.10  °C 
28.32  °C 
100.65  °C 
21.7% 


1 .  Determination  of  Average  Bulk  Temperature 


T 


b 


(T  .  +  T 
ci 


co 


)/2 


T.  =  (24.10  +  28 . 32)  /2 

D 


T,  =  26.21  °C 

D 


2.  Thermophysical  Properties  (evaluated  at  T^) 
Pr  =  5.83 

p„  =  997  kg/m3 
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855  x  10  N •  s/m 


Cpc  =  4.179  kJ/kg-K 


kc  =  618  xio-3  W/m*K 


mc  =  0.243  kg/s 


Determination  of  Cooling  Water  Velocity 


V  =  xn  /(p  *A. ) 
c  c'  i 


V  =  (0 . 243)/ (997) • (1. 56  xlO-4) 


Vc  =  1.56  m/s 


Determination  of  Reynolds  Number 


Re  =  (pc  Vc  Di)/pc 


Re  =  (997)  (1.56)  (.014D/855  xio 


-6 


Re  =  25,649 


6.  Determination  of  Heat  Transfer 


7.  Determination  of  Heat  Flux 


q  =  Q/U  Dq  L) 
q  =  4285/tt(0. 015875)  (0.305) 

q  =  281,700  W/m2 


8.  Determination  of  Nusselt  Coefficient 


0.651 


Of  hfQ  (9.81) 
Of  Do  * 


1/3 


Assume: 


T-  =  T  . 
f  sat 


Tf  =  100.65  °C 


=  957.8  kg/m' 


h  =  2255  J/kg 

kf  =  680.52  W/m-K 

Mf  =  277.6  xlO*6  N-s/m2 

„  .r-,  ,{680.52xl0“3)  3  (957.8)  2(2255xl03)  (9.81),1/3 

ftNn  =  - Tfi - J 

wu  (277.6x10  °)  (0.015875)  (281,700) 

hNu  =  11,243.6  W/m2 • K 
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9.  Determination  of  T, 


Tf,c  -  Tsat  -  <’/<2'hNu) 


T.  =  100.65  -  281,700/2(11,243.6) 

f  ,c 


T.  =  88.12  °C 

f  /  c 


10.  Thermophysical  Properties 


=  674.6  *  10  W/m*K 


=  966  kg/m' 


320.3  xio-6  N-s/m2 


=  2288.3  *  10°  J/kg 


11.  Determination  of  Nusselt  Coefficient 


"k 3  p2  h  9.81  ]  1/3 

V  *  O'651  - 


„  r(674.6xl0“3)3(966)2(2288.3xl03) (9.81),1/3 

U  .  03  i  l - 2 - J 

(320.3xl0_O) (0.015875) (281,700) 


h„  =  10,739.6  W/m  K 

Nu 


12.  Determination  of  T. 


Tf,o  *  Tsat  -  ^(2-hNu> 
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T _  =  100.65  -  281,700/(2) (10,739.6) 

£  9  C 

T,  =  87.53  °C 

f /C 


Determination  of  Logarithmic  Mean  Temperature 

T  -T  . 

LMTD  =  (Tcq  -  Tei)/in(^t--T^) 

sat  co 

LMTD  =  (28.32  -a4.1Q)/tn(^g;||4||.^) 

LMTD  =  74.42  °C 

Determination  of  Overall  Heat-Transfer  Coefficient 
UQ  =  q/LMTD 

UQ  =  281,700/74.42 

UQ  =  3785.3  W/m2 • K 

Determination  of  Inside  Heat-Transfer  Coefficient 


hi  =  ^|^|t(0.028)  (25, 649)°*8(5.83)°-333  (l.l) 

hL  =  8181  W/m2K 

Determination  of  inner  Wall  Temperature 

Tw  =  Tb  +  9  Do/(hi  Di) 

Tw  =  26.21  +  (281,700)  (0.015875)/ (8181)  (0.0141) 

T  =  64.98  °C 

w 


Determination  of  y  at  the  Average  Wall  Temperature 


y  =  433  xlO"6  N-s/m2 

w 


Determination  of  Correction  Factor 


C,.  =  (y  /y  ) 

f ,c  c'  w 


0.14 


Cf,c  =  1‘1 


Determination  of  Outisde  Heat-Transfer  Coefficient 


1U 

_ _  2  —  R 

U  D.  h.  w 
O  11 


I  0.015875 

3785.3  “  (0.0141)  (8181) 


-  0.000042925 


11,957  W/m2K 


APPENDIX  B 


UNCERTAINTY  ANALYSIS 

The  uncertainty  analysis  performed  for  this  thesis  is 
based  on  the  development  described  in  Kline  and  McClintock 
[Ref.  33] .  Kanakis  [Ref.  12]  employed  this  development  and 
derived  the  expressions  for  uncertainty  in  the  variables  of 
this  thesis  also.  In  this  thesis,  the  results  presented  by 
Kanakis  [Ref.  12]  in  Appendix  C  will  be  utilized,  but  where 
necessary,  corrections  will  be  pointed  out. 


A.  UNCERTAINTY  IN  THE  COOLING  WATER  VELOCITY 


6  V 
_ c 

V 


6mc  2  6pc  2  6Ai  2 

(-—r  +  (-—)  +  (-r^) 

m  pc  Ai 

c 


1/2 


with  the  following  uncertainties  assigned: 


Smc  =  ±0.01  kg/s 

5pc  =  ±3  kg/m3 

6A .  =  ±  0.00000022  m2 

l 


B.  UNCERTAINTY  IN  THE  REYNOLDS  NUMBER 


6  Re 
Re 


6pc  2  6Vc  2  5Di  2  6yc  2 
PC  Vc  i  Mc 


1/2 
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1  »  ■  *  ■  W  '*  I  1  f  » 


with  the  following  uncertainties  assigned: 


6pc  =  ±3  kg/m' 


SVc  as  calculated  in  Section  A 


W 

'-r. 


5Di  =  ±  0.0001  m 


=  ±  10  x 10  ®  N*s/m^ 


C.  UNCERTAINTY  IN  HEAT  TRANSFER 


6Q  _ 
Q 


6m_  ,  6  T  0  6  T  .  _  6C  0 

(~ ~)  2  +  (f~-T~  )2  +  <T— T  >  +  (-c^) 

mc  co  ci  co  ci  pc 


1/2 


with  the  following  uncertainties  assigned: 


6mc  =  ±0.01  kg/s 


6T  „  =  ±  0.025  °C 

co 


<5T  .  =  ±  0.025  °C 

ci 


&C  -  ±  2  J/kg-C 

pc 


D.  UNCERTAINTY  IN  THE  HEAT  FLUX 


£st  = 
q 


<^>2  -  (^>2  +  <^>2 


1/2 


with  the  following  uncertainties  assigned: 
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<$Q  as  calculated  in  Section  C 
6DQ  =  ±  0.0001  m 

6L  =  ±  0.001  m 
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F.  UNCERTAINTY  IN  T 


'f /C  _ 


5T  .  ?  o  6h»T  9  1  1/2 

(Hp— V  +  (^)2  +  (fr^)2 
sat  q  nNu 


with  the  following  uncertainties  assigned: 


5T  =  ±  0.025  °C 

sat 


6q  as  calculated  in  Section  D 


<5h„  as  calculated  in  Section  E 
Nu 


G.  UNCERTAINTY  IN  LOGARITHMIC-MEAN-TEMPERATURE  DIFFERENCE 


6LMTD 

LMTD 


{Tsat(Tco-Tci» 


^Tsat-Tci  ^  ^sat^co^  ^ 


T  -T  .  I 
sat  ci.  1 

T  -T  ’  / 


sat  co 


*( 


T  -T  . 

(T  -T  )  £n  (  saT: - — ) 

Usat  W^T  -T  1 
sat  co 


T  .  ~T  . 

(T  -T  )  Jin  (  sa*~ _ £i) 

Usat  ci  T  . -T 

sat  co 


with  the  following  uncertainties  assigned; 


)1 


l.*.  wV« -•-Ju’Lifc-V*-*** 


±  0.025  °C 


6T  =  ±  0.025  °C 

co 


6T  .  =  ±  0.025  °C 

ci 


H.  UNCERTAINTY  IN  OVERALL  HEAT-TRANSFER  COEFFICIENT 


,6q,2  .  ,SLMTD,2 
(  q  1  (  LMTD ’ 


with  the  following  uncertainties  assigned: 


<5q  as  calculated  in  Section  D 


6 LMTD  as  calculated  in  Section  G 


I.  UNCERTAINTY  IN  INSIDE  HEAT-TRANSFER  COEFFICIENT 


,  i.2  ,0.8  6  Re.  2  ,0.333  6Pr  .  2 

'  V  '  +  '  n  _  '  +  '  pg  '  +  ' 


6C.  y  0.146 (u/u)  2  1/2 

+  ( _ i.)  £  +  ( _ c  . 

+  '  C.  1  y  /y  ' 

l  Mc'  Hw  J 


with  the  following  uncertainties  assigned: 


6k  =  ±  0.0005  W/m'K 

c 


6D .  =  ±  0.0001  m 


6 Re  as  calculated  in  Section  B 


<5P  =  ±  0.063 

r 

6Ci  =  ±  0.0001 


6 (u  /y  )  =  ±  4.9  xlO-6  N-s/m2 

c  w 


J.  UNCERTAINTY  IN  TEMPERATURE  DIFFERENCE 


SDT 

DT 


.  0  5h.  -  0 

+  (i—)2  +  (-5s)2  +  ( 


with  the  following  uncertainties  assigned: 


Sq  as  calculated  in  Section  D 


Sh.  as  calculated  in  Section  I 
i 

6D  =  ±  0.0001  m 

o 

6D. 


±  0.0001  m 


For  the  data  presented  in  Appendix  A,  the  following  uncertain 
ties  are  the  result: 

quantity  uncertainty 

V  *  1.56  ±  0.064  m/s 


Re 


25,649  ±  1112 


APPENDIX  C 


DATA  REDUCTION  AND  PLOTTING  PROGRAMS 


1000!  FILE  NAME:  DRP3 
1005!  DATE:  Hay  9.  1984 

1010! 

1015  BEEF- 

1020  PRINTER  IS  1 

1025  PRINT  USING  "4X  .""Select  Option: . 

1030  PRINT  USING  "6X."”0  Taking  data  or  re-processing  previous:  data . 

1035  PRINT  USINC  "6X.""1  Plotting  previous  data . 

1040  PRINT  USING  "6X.'"’2  Computing  exponent  for  experimental  data . 

1045  PRINT  USING  "6X,""3  Plotting  on  LOG-LOG . 

1050  PRINT  USING  "6X,""4  Labelling . 

1055  INPUT  lop 

1060  IF  1op-0  THEN  CALL  Ham 

1065  IF  lop* I  THEN  CALL  Plot 

1070  IF  Iop*2  THEN  CALL  Expo 

1075  IF  lop*3  THEN  CALL  Lplot 

1080  IF  Iop-4  THEN  CALL  Label 

1 085  END 

1090  SUB  Hain 

1095  COM  /Ci/  C<7> 

1100  DIM  Tci t2) ,Tco<4,2> . T i <41 ,Mftt4) . Vu<  4> ,Ho<4) 

1105  DIM  To<4> . Ts<  1  > ,Tb<4) ,R3<4> .R4(4> ,S3<4> .S4<4> ,Fc(4» ,Mf i (4) .Mfs(4> 
1110! 

1115!  ASSIGN  COEFFICIENTS  FOR  THE  8-TH  ORDER 
1120!  POLYNOMIAL  FOR  TYPE-T  < COPPER -CONST ANT AN > 

1125!  THERMOCOUPLES 

1130  DATA  0.10086091 .25727.94369.-767345.8295.78025595.81 
1135  DATA  -9247486589 , 6 . 97688E ♦11.-2.661 92E +13,3. 940 78E+ 1 4 
1140  READ  C<») 

1145! 

H50!  ASSIGN  CONSTANTS  FOR  ROTAMETER  CALIBRATION 
1155!  LINES  (kg/miri) 

1160! 

1165  DATA  -1 -256E-4.5.866E-3.2.456E-4.-3.352E-3.-5.878E-4 
1170  DATA  1 . 121E-2 . 1 .237E-2, 1 .221E-2 . 1 .233E-2. 1 .232E-2 
1175  READ  Mfid.NfsC) 

1180* 

1185!  ASSIGN  SIEDER-TATE  COEFFICIENT  AND  EXPONENT 
1190!  FOR  REYNOLDS  NUMBER 
1 ’95  Ex*. 8 
1200! 

1205!  ASSIGN  GEOMETRIC  VARIABLES 
1210  Di*.0141  !  Inner  diameter  <m> 

1215  Do*. 015875  !  Outer  diameter  <m> 

1220  Ktm*21 .9  !  Thermal  conductivity  of  titanium  <W/r:-K> 

1225  L*.305  !  Condensing  length  <m> 

1230  Pt*Do*1 .5  !  Transverse  tube  pi tch-to-dia»eter  ratio 

1235! 

1240!  COMPUTE  THE  MEAN  STEAM  FLOW  AREA  IN  THE  TEST  CONDENSER  <m  2> 

1245  Af * (Pt ' 2-PI*Do  2/4)/Pt»2*L 
1 250 ! 

1255!  COMPUTE  INSIDE  AREA  AND  WALL  RESISTANCE 

1260  Ai*PI*Dr2/4 

1265  Ru*Do«L0G< Do/Di > /< 2»Ktm > 

1270! 

1275  PRINTER  IS  70' 

1280  CLEAR  70S 
1285  BEEP 

1290  INPUT  "ENTER  MONTH.  DATE.  AND  TIME  <MM : DD :HH :MM:SS>" . T imeS 
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1295  OUTPUT  709;"TD":Ti»«l 

1300  Sgy-0 

1305  Sgf-0 

1310  Sgyf-0 

1315  Sgf 2*0 

1320  BEEP 

1325  INPUT  "VAPOR  SHEAR  DATA  < 1 -Y . 0-N>?" , Iwd 
1330  IF  lvd-0  THEN 
1335  BEEP 

1340  INPUT  "ENTER  OUTPUT  NODE  <1 -SHORT .O-LONG)" . Jop 
1345  END  IF 
1350  BEEP 

1355  INPUT  "ENTER  THE  INPUT  NODE  ( 1-3054A. 2-FILE)". Im 
1360  IF  Im*2  THEN 
1365  BEEP 

1370  INPUT  "ENTER  THE  NAME  OF  THE  EXISTING  DATA  FILE",01ddataS 

1375  PRINT  USING  "10X.“"This  analysis  is  For  data  file  "" , 1 OA" :01ddatal 

1 380  BEEP 

1385  INPUT  "ENTER  1  IF  PROCESSING  1983  DATA" . Idt 

1390  ASSIGN  8File2  TO  OlddataS 

1395  END  IF 

1400  IF  In-1  THEN 

1405  BEEP 

1410  INPUT  "GIVE  A  NAME  FOR  THE  DATA  FILE  TO  BE  CREATED" .NewdataS 
1415  CREATE  BOAT  N  itaS  '*5 

1420  ASSIGN  #F  1 1  e '  Nei  .taS 

1425  END  IF 
1 430  BEEP 

1435  INPUT  "GIVE  A  NAME  FOR  THE  OUTPUT  FILE" .F i le_outS 
1440  BEEP 

1445  INPUT  "ENTER  TUBE  TYPE  (0-PLAIN. 1 -ROPED)" . I tt 
M50  IF  Itt-0  THEN 

1455  PRINT  USING  "10X.""Tube  type  :  Plain  . 

1460  Ci-. 026 

1465  END  IF 

1470  IF  Itt-1  THEN 

1475  PRINT  USING  "10X.""Tube  type  :  Roped  . 

1480  Ci*. 057 
1485  END  IF 

1490  PRINT  USING  " 1 0X . ""S ieder-Tate  constant  -  "".Z.4D":Ci 
1495  BEEP 

1500  INPUT  "ENTER  THE  INUNDATION  CONDITION  (1-5  TUBES.  2-30  TUBES)". Mi 
1505  IF  Mi-2  THEN  PRINT  ”  Inundation  condition  :  30  TUBES" 

1510  IF  Mi-1  THEN  PRINT  "  Inundation  condition  :  5  TUBES" 

1515  IF  Ivd-0  THEN 
1520  BEEP 

1525  INPUT  •■•WANT  TO  IN<  VAPOR  SHEAR  (  i  -  Y .  D-N)?"  .  Ivs 
1530  IF  I vs  *  1  THEN 
1535  BEEP 

1540  INPUT  "ENTER  B-VALUE  FOR  VAPOR  SHEAR  CORRELATION" . B 
1545  BEEP 

1550  INPUT  "ENTER  EX'  1NENT  FOR  VAPOR-SHEAR  CORRELATION"  .Nvs 
1555  END  IP 
1560  END  IF 

1565  CREATE  BDAT  rile_out$.6 

1570  ASSIGN  «File3  TO  File_outS- 

1575  Ja-0 

1580  Nrun*(l 

1585  FOR  1*0  TO  4 

1590  S3( I)-0. 

1595  S4( I >-0. 


1600  NEXT  I 
1605  IF  Im-1  THEN 
1610  BEEP 

1615  INPUT  "ENTER  MANOMETER  READINGS  <  HL  H . HL  H . HRU , HRH  > " . H 1 u . H 1 m , H  r  w . H  r  m 
1620  OUTPUT  •File1:Hlw,Hlm.Hrw,Hrm 
1 625  ELSE 

1630  IF  Idt-0  THEN  ENTER  9File2.-Hlw.Hlm.Hru.Hrm 
1635  END  IF 

1640  Dp*Hlm-Hrm+<Hlw-Hlm-Hru+Hrm>/13 .5 
1645  Mdot'FNMdot(Dp) 

1650  Nrun*Nrun+1 
1655  OUTPUT  709:"TD" 

1660  ENTER  709:TimeS 

1665  IF  Ivd-0  AND  Jop-0  THEN  PRINT  ••  •• 

1670  IF  Nrun-1  THEN 

1675  PRINT  USING  "10X .""Month,  date,  and  time:  "". 15A":Tn»e« 

1680  PRINT 
1 685  END  IF 
1690  IF  Im-2  THEN  Rdf 
1695  BEEP 

1700  INPUT  "ENTER  FLOW  METER  READINGS  (AS  PERCENTAGES)" .Fml .Fm2 
1705  DISP  “START  COLLECTING  CONDENSATE" 

1710  BEEP 
1715  WAIT  20 

1720  OUTPUT  709:"AR  AF0  AL19" 

1725  OUTPUT  722: "FI  R1  T1  Zl  FL1“ 

1730! 

1735!  READ  INLET  WATER  TEMPERATURES 
1740! 

1745  FOR  1-0  TO  2 
1750  OUTPUT  709: "AS  SA" 

1755  ENTER  722:Tci<I> 

1760  CALL  Tvsv( Tci < I ) ) 

1765  Tci < I )-FNTemp< Tc i < I ) . I ) 

1770  NEXT  I 
1  775* 

1780!  READ  OUTLET  WATER  TEMPERATURES 
1785! 

1790  I i-2 

1795  FOR  1*0  TO  * 

1800  IF  1-0  OR  1*3  THEN 

1805  Iu*2 

1810  ELSE 

1615  Iu-1 

1820  END  IF 

1825  FOR  J-0  TO  Iu 

1830  I  i  - 1 1  +  1 

1835  OUTPUT  709: "AS  SA" 

18*0  ENTER  722  :  Tco< I . J> 

1845  CALL  Twsv< Tco< I , J> ) 

1850  Tco< I , J) *FNTemp< Tco(I.J).Ii) 

1855  NEXT  j 
1860  NEXT  I 
1865! 

1970!  READ  STEAM  TEMPERATURES 
1875! 

1880  FOR  1-15  TO  16 
1885  OUTPUT  709: "AS  SA" 

1390  ENTER  722: Ts< 1-15* 

1895  CALL  Tvsv(Ts<T-15> > 

1900  Tf.(  I  -  15>-FNT emp<  Ts*  I  -  1 5 )  .  I  > 


1905  NEXT  I 
1910! 

1915!  READ  CONDENSATE  TEMPERATURE 
1920! 

1925  OUTPUT  709: "AS  SA" 

1930  ENTER  722:Tcon 
1935  CALL  TwsviTcon) 

1990  Tcon-FNTenpt Tcon , 1 7 ) 

1995! 

1950!  READ  VAPOR  TEMPERATURE 
1955! 

I960  OUTPUT  709; "AS  SA" 

1965  ENTER  722: Tv 
1970  CALL  TvsvUv) 

1975  Tv-FNTei»p<Tv.18> 

1980! 

1985!  READ  VAPOR  PRESSURE 
1990! 

1995  OUTPUT  709: "AS  SA" 

2000  ENTER  722:P_volts 
2005! 

2010!  COMPUTE  AVERAGE  HATER  TEMPERATURES  AT  INLET 
2015! 

2020  Ti<0)*Tci<0) 

2025  Ti<1 >-<Tci(0>+Tcj<1 >>-.5 
2030  T i<2)-Tci<  1  > 

2035  Ti<3)-<Tci<1)+Tci<2>>».5 

2090  Ti<9>*TCi<2> 

2095! 

2050!  COMPUTE  AVERAGE  WATER  TEMPERATURES  AT  OUTLET 
2055! 

2060  FOR  1-0  TO  9 
2065  IF  1-0  OR  1-3  THEN 

2070  To<I)-<Tco<I.0>-*Tco<I.1)  +  Tco<I.2>>».  3333 

2075  ELSE 

2080  Tq<I)-tTco<I.0)+TcotI.1>)«.5 

2085  END  IF 
2090  NEXT  I 

2095  Tsa-<Ts<0>+Ts< 1 ) )« .5 
2100  Pvap-FNPvsv(P_vol ts> 

2505  Tsat*FNTvsp<Pvap»1 33 . 322 > 

2110  Dsup*Tv-Tsat 
2115! 

2' 20!  READ  INFORMATION  FOR  CONDENSATE  FLOW  RATE 
2125! 

2130  IT  Ivd-0  THEN 
2135  BEEP 

2590  INPUT  "ENTER  INITIAL  AND  FINAL  LEVELS  IN  HOT  WELL  1".H1.H2 
2195  DH-H2-H1 

2 '50  IF  Nrur.  MOD  5-1  THEN  Nsum-O 
2155  Mf 1 *590 .9836»Dh 
2 : 60  Mdl-Mf 1»FNRhow<Tsat-10)«1 .0E-6/60 
2165  Msujn*Msum+Mf  1 

2' 70  IF  Mi -2  AND  Nrun<>30  AND  Nrur,  MOD  5-0  THEN 
2175  Mave-Mstin/5 

2180  Set*  <Mave«FNRhow<  Tsat-  10) / 1 0  6-*  .03238) /. 0921  32 

2185  END  IF 

2190! 

2195  Rdf:  ! 

2200! 

2205  IF  lvd*0  AND  Jop*0  THEN 


2210 
2215 
2220 
2225 
2230 
2235 
2240 
2245 
2250 
2255 
2260 
2265 
2270 
2275 
2280 
2285 
2290 
2295 
2300 
2305! 

2310!  CALCULATE  AVERAGE  BULK  TEMPERATURES 
2315! 

2320  IF  Ivd-0  THEN  Nx-4 

2325  IF  Ivd-1  THEN  Nx-0 

2330  FOR  1-0  TO  Nx 

2335  Tb<I>*(Ti(I)+To<I>>«.5 

2340  NEXT  I 

2345! 

2350  IF  Mi  *  1  OR  (Mi-2  AND  Nrun<6>  THEN  Asl-0. 

2355  IF  Mi-2  AND  Nrun>S  THEN  S'-Asl 

2360  Sl-Asl 

2365  FOR  J-0  TO  Nx 

2370  IF  J-0  THEN  Cwf-Fml 

2375  IF  J-1  THEN  Cuf-Fm2 

2380  Mf-Mf i(J>+Mfs(J>*Cuf 

2385  Tx*Tb( J  > 

2390  Vu<J>-Mf/(FNRhow(Tx>»Ai) 

2395' 

2400!  CALCULATE  INSIDE  AND  OUTSIDE  COEFFICIENTS 
2405! 

2410  Rew-FNRhou(  T  x  >»Vu<  J)«Di /FNMuw(  T x ) 

2415  Cf-1. 

2420  Q-Mf  «FNCpu<  Tx  >»(To(J>-Ti<J)> 

2425  IF  Nrun-1  AND  J-0  THEN  01 -Q 
24 30  Qp»Q/(PI*Do«L ) 

2435  IF  <Mi- 1  OR  <Mi-2  AND  Nrun<6)>  AND  J-0  THEN 
2440  T  f i lm-Tsat 
2445  K.f -FNKut  Tf  i  lm) 

2450  Rhof -FNRhou(  Tf i lm  > 

2455  Hfg»FNHfg(Tsat>*1000 
2460  Muf-FNMuw<  Tf  1 1  m  > 

2465  Hnu-  .651 -Kf  *< Rhof  '  2-Hf g*9 . 799/ (Muf -Do-Op >  >  .3333 

2470  Tf i lmc-Tsat-Qp/Hnu* .5 

2475  IF  A6$<  <  Tf i Imc- Tf ilm)/Tf i Inc  >  > . 01  THEN 

2480  Tfilm-Tfiimc 

2485  GOTO  2445 

2490  END  IF 

2495  IF  Ivd-0  AND  Jop-0  THEN 

2500  PRINT  USING  "10X,""Nusselt  coefficient  for  first  tube  -  '"'.5D.D.'"'  <W/m  ? 

K> . :Hnu 

2505  END  IF 


PRINT  USING  "10X,""Run  number  *  "".DD“:Nrun 

PRINT  "  Tube  *  1  2  3  4  5 

END  IF 
END  IF 

IF  Ivd-1  AND  Nrun-1  THEN 

PRINT  USING  "1  OX. ""Data  •  Vv  Re-tp  Nuc  Y  F . 

END  IF 

IF  Im-2  THEN 

IF  Nrun  MOD  5-1  AND  Mi-2  AND  Nrun>5  THEN  ENTER  «File2:Fpt 
ENTER  #File2:Ti<»),To<*>.Tsa.Tcon. Tv .Pvap .Tsat . Dsup . Fn 1 ,Fm2 
ENTER  8File2:H1 ,H2 
END  IF 

IF  Ivd-0  AND  Jop-0  THEN 

PRINT  USING  "10X. ""Inlet  temp  <Deg  C>  :"".5(DDD.DD.2X)";Ti(«> 

PRINT  USING  *•  1 0X .""Out  let  temp  (Deg  C> :"".5(DDD.DD.2X>";To(-> 

PRINT  USING  "10X. ""Saturation  temperature  -  "".3D.DD."”  (Deg  C> . ; Tsat 

PRINT  USING  "10X. “"Degree  of  superheat  -  "".3D.DD.”"  (Deg  C> . ;Dsup 

PRINT  USING  ”10X. ""Static  pressure  -  "“.3D.DD.""  (mm  Hg> . :Pvap 

END  IF 
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2510  END  IF 

2515  If  J-0  AND  Ivd-0  AND  Jop-0  THEN 
2520  IF  Mi  - 1  AND  Nruu-1  THEN  Ho1*0. 

2525  PRINT  “  Tube  Vw  Heat  f lu.*  Cono  coef  R1 

RR" 

2530  PRINT  *’  •  <m/s>  ‘U/m  2>  (W/m  2.K>" 

2535  END  IF 

2540  If  Nruri*  t  AND  Jvs-1  AND  Ivd-0  THEN 
2545  Vg-FNVvsttTsat) 

2550  Vv-Mdot-Vg/Af 
2555  Dt*<Tsat-Tf iln>«2 
2560  Muf-FNMuw<Tf  ill*) 

2565  Hfg-FNHfg(Tsat >-1000 
2570  Kf*FNKw<Tf ilm) 

2575  Rhof-FNRhow<Tf  1  li*> 

2580  Retp-Rhof -Vv-Do/Huf 

2585  F  f -9 . 799-Do-Muf  *H  f  g  /  ( Vv '  2  -K  f  -D  t  > 

2590  IF  Ivd-0  THEN 

2595  Nuc-B-Ff  Nvs-Retp  .5 

2600  Hoc-Nuc-Kf/Do 

2605  Tf  i  Iinc-Tsat-Qp/Hoc- .5 

2610  IF  AB$<Tf iIi*-Tfilmc)>.1  THEN 

2615  TF  j lm*Tf i Inc 

2620  GOTO  2555 

2625  END  IF 

2630  END  IF 

2635  END  IF 

2640  MuwFNMuul  T* ) 

2645  Hi  -FNKw< T *  > /Di -C i -Reu  Ex- (F NPrwf Tx ) >  * . 3333-CF 

2650  Dt*Qp7Hi-Do/Di 

2655  Cf  c- <Muu7<FNHuu< Tx+Dt>>)' . 14 

2660  IF  ABS<<Cf-Cfc>/Cfc>>.01  THEN 

2665  Cf-<Cf *Cfc>*.5 

2670  GOTO  2645 

2675  END  IF 

2680  Lmtd*<To< J>-Tit J)>/L0G((Tsat-Ti < J> ) /( Tsat-Tot J>>> 

2685  IF  Nrun-1  AND  Ivs-1  AND  Ivd-0  THEN 

2690  Fc<J>-Hoc/Hr.u*<0/Ql )  .3333 

2695  Mdot-Hdot-Q/Hf  g 

2700  END  IF 

2705  Uo-Op/Lntd 

2710  Ho<  J>  *  1 /<  1 /Uo-Do/<Di -hi ) -Ru> 

2715  IF  Ivd-1  THEN 
2720  Dt*Qp/Ho<  0  > 

2725  Tf ilm-Tsat-Dt/2 
2730  Og*FNVvst < Tsat > 

2735  Rhof-FNRhow<Tf  ili*> 

2740  Vv-Mdot-Vg/Af 
2745  Muf-FNIWTf  ilm) 

2750  Kf *FNKw< Tf i lm> 

2755  Hfg-FNHf  g<  Tsat  >-1000 

2760  Retp*Rhof -Vv-Do/Huf 

2765  Ff-9.799-0o-Muf*Hfg/<Vv  2-Kf-Dt> 

2770  Nuc-Ho(0)-Do/Kf 
2775  Y-lOG<Nuc/Retp  .5) 

2?80  F-L0G<Ff) 

2785  Sgy-Sgy4Y 
2790  $gf-$gf  *F 
2 795  Sgyf -Sgyf+T-F 
2900  Sof2*Sp1?4F  2 
2805  END  IF* 
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2b1U  if  lvd-0  THIN 

2815  IF  lv«-1  THEN  Ho< J)-Ho< J)/Fc< J) 

2820  Rr*Uo/Ho< J) 

2825  S1-S1+Ho<J> 

2830  IF  Nrun  MOD  5-1  THEN 

2835  IF  Ni-1  OR  <Mi-2  AND  Nrun-31>  THEN  Ja-0 

2840  IF  Mi -2  AND  5<Nrun  HND  Nran<30  THEN  Ja-Nrun-1 

2845  IF  Mi-2  AND  35<Nrun  THEN  Ja-Nrun-1 

2850  END  IF 

2855  IF  Mi-1  OR  <30<Nrun  AND  Nrun<36  AND  Mi  - 1 >  OR  Nrun<6  THEN 
2860  R1-Ho<J)/Ho<0> 

2865  R2-S1/<<J+1+Ja>«Ho(0>> 

2870  ELSE 

2875  R1*Ho(J)/Ho1 

2880  R2-S1 /<  <  J+ 1 + Ja>»Ho1 ) 

2885  END  IF 
2890! 

2895!  PRINT  RESULTS 
2900! 

2905  IF  Jop-0  THEN 

2910  PRINT  USING  "1 1 X .DD .4X ,DD .DD ,2X ,2<D .5DE ,?X> ,3<Z .4D.2X)" : J+ t+Ja.Vu( J> .Qp.Ho 
<J).R1 .R2.Rr 

2915  IF  Im-1  AND  J-4  THEN 
2920  BEEP 

2925  INPUT  "OK  TO  ACCEPT  THIS  DATA  SET  <  1  -Y  ,0-N>'?”  .Oke 

2930  IF  Ok  sOI  THEN 

2935  Nrun-Nrun-1 

2940  GOTO  1650 

294S  END  IF 

2950  END  IF 

2955  END  IF 

2960! 

2965  IF  Mi -2  AND  Nrun<6  AND  J-0  THEN 

2970  Hoi  -Hoi +Ho<  0  > /5 

2975  END  IF 

2980  FOR  K-0  TO  4 

2985  IF  K-J  THEN  S3<K)-S3(K>+R1 

2990  IF  K-J  THEN  S4<K)-S4<K i+R2 

2995  NEXT  K 

3000  END  IF 

3005  NEXT  J 

3010  IF  Ivd-0  THEN 

3015  IF  Nrun  MOD  5-0  THEN 

3020  FOR  K-0  TO  4 

3025  R3(K)-S3<K)/5 

3030  R4<K)-S4(K)/5 

3035  S3<K>-0. 

3040  S4(K)-C. 

3045  NEXT  K 

3050  IF  Hi-2  AND  Nrun  MOD  5-0  AND  Nrurt<>3D  THEN  Asl -Nrun«R4(4>«H01 
3055* 

3060!  PRINT  AVERAGE  RATIOS 
3065! 

3070  IF  Jop»0  OR  <Jop*1  AND  Nrun<6>  THEN 
3075  PRINT 

3080  PRINT  "  Tube  *  R3  R4" 

3085  END  IF 
3 uSO  TOR  J-i  TO  5 

3095  PRINT  USING  "12X . DD,2<4X . Z .4D>" : J  + Ja.R3< J- 1 > ,R4<  J- 1  > 

3100  OUTPUT  6FiIe3:J-*Ja.R3<J-l  >  ,R4(J-1  ) 

1 05  NEXT  J 
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3110  END  IF 

3115  IF  Nrun  HDD  5*0  AND  Mi-2  AND  Nrun<>30  AND  Im-1  THEN 
3120  BEEP 

3125  PRINT  USING  "10X.""Set  porous-tube  flowmeter  reading  to  "".3D.D.""  PERCENT 

. :Set 

3130  END  IF 
2135  END  IF 
3140  IF  Ivd-1  THEN 
3145  Ey-EXP(Y) 

3150  PRINT  USING  "1 2X .DD .5X , Z . DD .2X . Z . 3DE .2X . 3D . DD .2X ,2< Z . 3D .2X >" ;Nrun , Vv .Retp . 
Nuc  ,Ey  .Ff 

3155  IF  Im-1  THEN 
3160  BEEP 

3165  INPUT  "OK  TO  ACCEPT  THIS  DATA  SET  d - Y . 0-N> 7" .Ok s 

3170  IF  OksOI  THEN 

3175  Nrun-Nrun* 1 

3180  GOTO  1605 

3185  END  IF 

3190  END  IF 

3195  OUTPUT  SFile3;Ff.Ey 

3200  END  IF 

32051 

3210  IF  Nrun  MOD  5-1  AND  Mi-2  AND  Nrun>5  AND  Im-1  THEN 

3215  OUTPUT  BFi lei :Set 

3220  Mpt--8. 3616! 3+ 10. 076742-Set 

3225  END  IF 

3230  IF  Im-1  THEN 

3235  OUTPUT  BF i lei ;Ti<->.To<-).Tsa.Tcon, Tv .  Pvap . Tsat ,  Dsup .Fml .Fm2 
3240  OUTPUT  SFilel :H1 .H2 
3245  BEEP 

3250  INPUT  "HILL  THERE  BE  ANOTHER  RUN  < 1  - Y . 0-N>7" ,Go_on 

3255  IF  Go_on- 1  THEN 

3260  IF  Ivd-0  THEN  1650 

3265  IF  Ivd-1  THEN  1605 

3270  END  IF 

3275  ELSE 

3280  IF  Mi-2  AND  Nrun<30  THEN 
3285  IF  Ivd-0  THEN  1650 
3290  IF  Ivd-1  THEN  1605 
3295  END  IF 

2300  IF  Mi  *  1  AND  Nrun<10  THEN 
3305  IF  Ivd-0  THEN  1650 
3310  IF  Ivd-1  THEN  1605 
3315  END  IF 
3320  END  IF 

3325  IF  Im-'  THEN  PRINT  USING  "10X.DD.""  Data  runs  were  stored  in  file  "".IDA": 
Nrun .NewdataS 
3330  PRINT 

3335  PRINT  USING  " 1  OX ,”"P 1 ot  data  are  stored  in  file  . 14A" : F i le  outs 
3340  IF  Ivd-1  THEN 

3345  Nvs* ( Nrun-Sgy f -Sgy-Sgf ) /< Nrun-Sgf 2-Sgf  2> 

3350  B-EXP< (Sgy-Nvs-Sgf )/Nrun> 

3355  PRINT 

3360  PRINT  USING  "10X. ""Vapor-Shear  Correlation: . 

3365  PRINT  USING  "12X.”"E  -  "”.Z.3D":f 

3370  PRINT  USING  "12X.""n  -  "" . Z . 3D" : Nvs 

3375  END  IF 

3380  ASSIGN  SFilel  TO  - 
3385  ASSIGN  «Fiie2  TO  * 

3390  ASSIGN  *File3  TO  . 


3395  SUBEND 
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3400! 

3405!  THIS  SUROUTINE  CONVERTES  THERMOCOUPLE  VOLTAGE  INTO  TEMPERATURE 
3410! 

3415  SUB  Tvsv<T> 

3420  COM  /C»/  C<7> 

3425  Suft’O . 

3430  FOR  1-0  TO  7 

3435  Stt»*Su»*C<I)-T  I 

3440  NEXT  I 

3445  T-Suft 

3450  SUBEND 

3455! 

3460!  THIS  FUNCTION  CALCULATES  PRANDTL  NUMBER  OF  HATER  IN  THE 

3465!  RANGE  15  TO  45  DEG  C 

3470! 

3475  DEF  FNPrw<T> 

3480  Y-10  <1 .09976605-TX1 .3749326E-2-TM3.968875E-5-3.45026E-7«T>>> 

3485  RETURN  Y 

3490  FNEND  • 

3495! 

3500!  THIS  FUNCTION  CALCULATES  THERMAL  CONDUCTIVITY  OF  HATER 

3505!  IN  THE  RANGE  OF  15  TO  105  DEG  C 

3510! 

3515  DEF  FNK*<T> 

3520  Y- .5625894+T»<2.2964546E-3-T«< 1 . 509766E -5-4 .0581 652E-8»T>  > 

3525  RETURN  Y 
3530  FNEND 
3535! 

3540!  THIS  FUNCTION  CALCULATES  SPECIFIC  HEAT  OF  HATER 

3545!  IN  THE  RANGE  15  TO  45  DEG  C 

3550! 

3555  DEF  FNCp«<T) 

3560  Y«<4.21 120858-T*<2 .26826E-3-T«<4 .42361E-5+2. 71428E-7«T>  >  >»1000 

3565  RETURN  Y 
3570  FNEND 
3575! 

35801  THIS  FUNCTION  CALCULATES  DENSITY  OF  HATER  IN  THE 

3585!  RANGE  15  TO  105  DEG  C 

3590! 

3595  DEF  FNRhow«T> 

3600  Ro-999.52946*T-< ,01269-T»<5.482513E-3-T*1 .234147E-5>> 

3605  RETURN  Ro 
3610  FNEND 
3615! 

3620!  THIS  FUNCTION  APPLIES  CORRECTIONS  TO  THERMOCOUPLE  READINGS 
3625! 

3630  DEF  FNTemp( T . I ) 

3635  DIM  A< 14)  ,B< 14) 

3640  DATA  0.640533.0.573054.0.593101 ,0.57298. 0.56228. 0.567384. 0.569577 

3645  DATA  0.553951 ,0.552008.0.566955.0.520998.0.522661 .0.531008.0.560788.0.5524 

05 

3650  DATA  1 1 . 8744 . 8 . 631 63 . 9 . 394 1 2 . 8 . 570246 . 8 . 299436 .8.36677.8.04507.7. 459766 
3655  DATA  7 . 498928 . 7 . 9408 . 5 . 870 72 . 5 . 391 556 . 6 . 1 3399 . 6 . 48586.6. 326224 
3660  READ  A<*).B(-> 

3665  IF  I < 1 5  THEN 

3670  T*T-<A<I)-B<I>-.001«T> 

3675  ELSE 
3680  T-T-.5 

3685  END  IF 
3690  RCTURN  T 
3695  FNEND 
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3700! 

3705!  THIS  FUNCTION  COMPUTES  THE  SPECIFIC  VOLUME  OF  STEAM 
3710! 

3715  DEF  FNVvst(Tt) 

3720  P-FNPvst<Tt> 

3725  T-Tt*273.15 

3730  X-1500/T 

3735  F1«1/<1*T«1.E-4> 

3740  F2-<1-EXP<-X>>'2.5*EXP<X>/X  .5 
3745  B- . 0015-F 1  -  . 000942-F2- . 0004882-X 
3750  K-2»P/<461.S2-T> 

3755  V-< 1+< 1 ♦2"B»K) ' .5>/K 
3760  RETURN  V 
3765  FNEND 
3770! 

3775!  THIS  FUNCTION  CONVERTS  THE  VOLTAGE  READING  OF  THE  PRESSURE 

3780!  TRANSDUCER  INTO  PRESSURE  IN  MM  HG 

3785! 

3790  DEF  FNPvsw< V)  . 

3795  Y-1 .1 103462063. 36413-V 

3800  RETURN  Y 

3805  FNEND 

3810! 

3815!  THIS  FUNCTION  CALCULATES  THE  SATURATION  TEMPERATURE  OF  STEAM  AS  A  FUNCTION 

3820!  OF  PRESSURE 

3825! 

3830  DEF  FNTvsp<P> 

3835  T  u« 1 1 0 
3840  T 1-80 
3845  Ta*<Tu+Tl >•  .5 
3850  Pc’FNPvst(Ta) 

3855  IF  ABS<<P-Pc>/P>> .0001  THEN 

3860  IF  Pc<P  THEN  Tl-Ta 

3865  IF  Pc>P  THEN  Tu-Ta 

3870  GOTO  3845 

3875  END  IF 

3880  RETURN  Ta 

3885  FNEND 

3890! 

3895!  THIS  FUNCTION  COMPUTES  THE  VISCOSITY  OF  HATER 
3900! 

3905  DFr  FNMuw<  T ) 

3910  A-247.8/<T*133.15> 

3915  Ma*2.4E-5«10'A 
3920  RETURN  Mu 
3925  FNEND 
3930! 

3935!  THIS  FUNCTION  COMPUTES  THE  LATENT  HEAT  OF  VAPORIZATION 
3940! 

3945  DEF  FNHf g<  T  > 

3950  Hfg-2497.7389-T«<r.2074+T-<l . 7079E- 3-2 . 8593E-6*T > > 

3955  RETURN  Hfg 
3960  FNEND 
3965* 

3970!  THIS  FUNCTION  COMPUTES  THE  SATURATION  PRESSURE 
3975! 


3980  DEF  FNPvst < Tstearn) 
3985  DIM  K<8> 

3990  DATA  -7.691234564. 
3995  DATA  4.16711732.20 


26 . 08023696 .-168.1 706546.64 
9750676. 1E9. 6 


23285504.-118.9646225 


4UUU  KtHb  M»  J 

4005  T-<Tstpa»«-273.15>/647.3 

4010  Su**0 

4015  FOR  N-0  TO  4 

4020  Su»*Su»+K<N)*<  1-T)‘<N-»1 ) 

4025  NEXT  N 

4030  Br-Su*/<T*<1+K<5>-<1-T>*K<6>-<1-T>‘2>>-<1-T>/<K<7>«n-T>~2+K<8>> 
4035  Pr-EXP(Br) 

4040  P-22120000-Pr 
4045  RETURN.  P 
4050  FNEND 
4055! 

4060!  THIS  FUNCTION  CALCULATES  THE  STEAM  MASS  FLOH  RATE 
4065! 

4070  DEF  FNMdot(Dp) 

4075  Mdot-4.3183E-3+5.6621E-4-Dp 
4080  RETURN  Mdot 
4085  FNEND 


APPENDIX  D 


SEIDER-TATE  CONSTANT  CALCULATION 

The  following  solution  procedure  used  in  calculating 
the  Sieder-Tate  constant,  by  using  a  modified  Wilson  plot 
is  presented  in  this  appendix: 


By  plotting  Z  vs.  W  and  computing  the  slope  of  the 
least-squares-fit  line,  and  then  taking  the  recip¬ 
rocal  of  this  slope,  the  resulting  value  is  the 
Sieder-Tate  constant. 


Tw  and  T£  must  be  determined  iteratively,  as 
described  in  Chapter  IV,  section  E. 
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